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Abstract 

We revisit an old idea that gravity can be unified with Yang-Mills theory by enlarging 
the gauge group of gravity formulated as gauge theory. Our starting point is an action that 
describes a generally covariant gauge theory for a group G. The Minkowski background 
breaks the gauge group by selecting in it a preferred gravitational SU(2) subgroup. We 
expand the action around this background and find the spectrum of linearized theory to 
consist of the usual gravitons plus Yang-Mills fields charged under the centralizer of the 
SU(2) in G. In addition, there is a set of Higgs fields that are charged both under the 
gravitational and Yang-Mills subgroups. These fields are generically massive and interact 
with both gravity and Yang-Mills sector in the standard way. The arising interaction of 
the Yang-Mills sector with gravity is also standard. Parameters such as the Yang-Mills 
coupling constant and Higgs mass arise from the potential function defining the theory. 
Both are realistic in the sense explained in the paper. 



1 Introduction 

There have been numerous attempts to unify Einstein's theory of gravity with gauge fields 
describing other interactions. One such unification proposal is that of Kaluza-Klein, where the 
metric and gauge fields arise from a higher- dimensional metric tensor upon compactification 
of extra dimensions. This scenario has become an indispensable part of string theory, which 
also provides another unifying perspective by viewing gravity and Yang-Mills as excitations of 
closed and open strings respectively. For more details on string-inspired unification schemes 
see a recent exposition [TJ. 
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There have also been attempts to unify gravity with gauge theory without introducing extra 
dimensions. There is, however, a very strong no-go theorem [2] that shows that at least one type 
of such unification is impossible. The theorem states that the symmetry group of the S-matrix 
of a consistent quantum field theory (in Minkowski spacetime) is the product of Poincare and 
internal gauge group. In other words, the spacetime and internal symmetries do not mix. The 
only way to go around this statement is via supersymmetric extensions of Poincare group [3]. 

Now, since gravity can be (at least loosely) viewed as a gauge theory for the diffeomor- 
phism group, and the later contains Poincare group as that of rigid global transformations, the 
Coleman- Mandula theorem [2] is sometimes interpreted as saying that no unification of gravity 
and gauge theory that puts together diffeomorphisms and gauge transformations is possible. In 
this discussion, however, one must be careful to distinguish between local gauge invariances of 
a theory and global symmetries whose presence or absence depends on a particular state one 
works with, see [I] that emphasizes this point. 

While it may be difficult or impossible to "unify" diffeomorphisms and gauge transforma- 
tions into a single gauge group, this is not the only possible way to approach the unification 
problem. To understand how a different type of unification might be possible, let us recall that 
in the so-called first-order formalism gravity becomes a theory of metrics as well as Lorentz 
group spin connections. The "internal" Lorentz group acts by rotating the frame and has no 
effect on the metric. Thus, the physical dynamical variable is still the metric, one simply added 
some gauge variables and enlarged the gauge group, which in this formulation is a (semi-) direct 
product of the diffeomorphism group and SO(l,3). Further, in the Hamiltonian formulation 
this theory can be easily cast into one on the Yang-Mills phase space. This is done by adding 
to the action a term that vanishes on-shell [5]. The phase space is then that of pairs SU(2) 
connection plus the canonically conjugate "electric" field. Thus, after the trick of adding an on- 
shell unimportant term, gravity becomes a generally covariant theory of an SU(2) connection. 
The spacetime metric (tetrad) is still a dynamical variable but in this formulation it receives 
the interpretation of the momentum canonically conjugate to the connection. 

Yang-Mills theory, on the other hand, after it is written for a general spacetime metric, also 
becomes a generally covariant theory of a connection and spacetime metric. One could then 
attempt to put the two generally covariant gauge theories together in some way that combines 
the "internal" gauge groups, while leaving the total gauge group to be a (semi-) direct product 
of diffeomorphisms and "internal" symmetries. This would not be in any conflict with the no-go 
theorem [2] for what is unified is not the Poincare and internal symmetry groups. This might 
not be what can be legitimately called a unification, for the end gauge group is not simple, but 
this idea does lead to some interesting "unified" theories, as we hope to be able to demonstrate 
in this paper. 

As far as we are aware, the first proposal of this type was put forward in [6J [7] , with the idea 
being precisely to extend the gauge group of gravity formulated in tetrad first-order formalism 
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as a theory of the Lorentz connection. This proposal was later pushed forward in [SJ [9], see 
also [10] for the most recent development. The key point of this proposal is that it is a non- 
degenerate metric that breaks the gauge symmetry of the unified theory down to a smaller 
group consisting of SO(l,3) for gravity and some "internal" group for Yang-Mills fields. 

A similar in spirit, but very different in the realization idea was proposed in [TT], and further 
developed in [121 EI3J EH]- This approach stems from the fact that Einstein's general relativity 
(GR) can be reformulated as a theory on the Yang-Mills phase space. At the time of writing 
[TT] it was achieved in Ashtekar's Hamiltonian formulation of GR [15] that interprets gravity as 
a special generally covariant (complexified) SU(2) gauge theory. The fact that gravity in this 
formulation becomes a theory of connection suggests that a gauge group larger than SU(2) can 
be considered. This is what was attempted in [HJ |T2l [131 E] > with the main result of [H] being 
that Yang-Mills theory arises in an expansion of the theory around the de Sitter background. 

The idea to put together the "internal" gauge groups of gravity and gauge theory is an 
interesting one. However, its particular realizations available in the literature are not without 
problems. Thus, the approach reviewed and further developed in [10] does a very good job in 
describing the fermionic content of the theory. Bosons, on the other hand, are described less 
convincingly in that many new propagating degrees of freedom (DOF) are introduced. The 
other approach |T4] is also not very convincing since it works at the phase space level, and it is 
generally very difficult to approach a theory if no action principle is prescribed. Another aspect 
of the particular realization given in [14] is that it naturally describes a complexified GR put 
together with complexified Yang-Mills. No natural reality conditions that would convert this 
into a physical theory were given. 

The unification by enlarging the internal gauge group proposal was recently revisited in 
[T6] . where the new action principle [TT] for a class of modified gravity theories [18] . extended 
to a larger gauge group was used. This work also avoided the reality conditions problem by 
extending the gauge group of an explicitly real formulation of gravity that works with the 
Lorentz, not with the complexified rotation group. Specifically, it was suggested in [16] that 
the action of the type proposed in [TT] considered for a general Lie group G describes gravity in 
its SO (4) part plus Yang-Mills fields in the remaining quotient G/SO(4). As in [TJ], the Yang- 
Mills coupling constant is related in [16] to the cosmological constant. As in the approach P,[T], 
in [16] it is a non-degenerate metric that breaks the symmetry down to a smaller gauge group. 
The approach of [16] is also similar to that of [61 [TJ in that many new bosonic degrees of freedom 
are introduced. Thus, it was shown in [19] that the BF-type action of [TT] for G = SO (4) does 
not anymore describe pure gravity theory in that it describes six new DOF. 

In this paper we take the described unification idea one step further. Our approach is 
similar in spirit to [16J in that we start from an action principle of the type first proposed in 
[TT] . However, unlike in [16], we interpret only a (complexified) SU(2) subgroup of the gauge 
group G as that corresponding to gravity. The part of the gauge group that commutes with 



3 



this gravitational SU(2) is then seen to describe Yang-Mills fields, and the part that does not 
commute with SU(2) describes charged scalar, i.e. Higgs, fields. We note that the suggestion 
that in unifications of this type the "off-diagonal" part of the Lie algebra corresponds to Higgs 
fields is contained already in [T6] . 

Our approach is also similar to the original proposal [H] that enlarged the SU(2) gravi- 
tational gauge group. However, in contrast to [H] that worked at the phase space level our 
starting point is an action principle that makes a much more systematic analysis possible. Also 
the details of our proposal differ significantly from that of [14J in that a semi-realistic (more 
on this below) unification is achieved without the need for a cosmo logical constant. Thus, the 
Yang-Mills coupling constant in our scheme is related not to the cosmo logical constant, which 
we set to zero, but to a certain other parameter of the theory. This features of our proposal 
also makes it different from that of [T6] . 

More specifically, we start from a generally covariant gauge theory for a (complex) semi- 
simple Lie group G, with certain reality conditions later imposed to select real physical con- 
figurations. A particularly simple solution of the theory describes Minkowski spacetime. This 
solution breaks G down to a (complexified) SU(2) times the centralizer of SU(2) in G. The 
spectrum of linearized theory around the Minkowski background is then shown to consist of the 
usual gravitons with their two propagating DOF, gauge bosons charged under the centralizer 
of SU(2) in G, and a set of scalar Higgs fields. The Higgs fields are in general massive, with the 
mass being related to a certain parameter of the potential defining the theory. After the reality 
conditions are imposed all sectors of the theory have a positive-definite Hamiltonian. We also 
work out interactions to cubic order and show that all interactions are precisely as expected. 
That is, all non-gravitational fields interact with gravity via their stress-energy tensor, and the 
interactions in the non-gravitational sector are also standard and are as expected for Higgs 
fields. Thus, our unification scheme passes the zeroth order test of being not in any obvious 
contradiction with observations. However, to obtain a truly realistic unification model many 
problems have to be solved. In particular, fermionic DOF are not considered in this paper 
at all. Thus, our results provide only one of the first steps along this potentially interesting 
research direction. We return to open questions of our approach in the discussion section. 

In this paper we have illustrated the general G case by considering the simplest non-trivial 
example of G = SU(3). This example is rather generic, and the same technology that we 
develop for G = SU(3) can be used for any Lie group. We could have presented a general 
semi-simple case treatment phrased in terms of the root basis in the Lie algebra. However, at 
this stage of the development of the theory it is not clear if there is any added value in doing 
things in full generality. We thus decided to keep our discussion as simple as possible and treat 
one example that, if necessary, is easily extendible to the general situation. 

Another general remark on this paper is as follows. As the reader will undoubtedly notice, 
a sizable part of our paper is occupied by the Hamiltonian analysis of various sectors, or of 



4 



the full theory. We also always give the Lagrangian treatment in which things are much more 
transparent. Thus, it might at first sight seem that the Hamiltonian formulation only clatters 
the exposition. We, however, believe that some aspects of the theory are much clearer precisely 
in the Hamiltonian formulation. For instance, our treatment of the reality conditions heavily 
uses the Hamiltonian analysis and it would be very hard to arrive at the correct conditions 
without it. This is our main reason for carrying out such an analysis in all cases that are 
discussed. 

The organization of the paper is as follows. In section [2] we define the class of generally- 
covariant gauge theories that is the subject of this paper. Section [3] performs a Legendre 
transformation that introduces the two-form field as the main dynamical variable and rewrites 
the action of our theory in a form most useful for practical computations. In section H] we 
sketch the Hamiltonian analysis and count the number of propagating DOF. Section [5] contains 
a general discussion on the problem of linearization. In section [6] we warm up by considering 
the case of pure gravity corresponding to G — SU(2). The Minkowski space background that 
we expand about is described here. Section [7] deals with an example of a non-trivial group for 
which we take G = SU(3). It is here that we obtain a Lagrangian describing the YM and Higgs 
sectors of our model. In section [8] we deduce interactions between various sectors of our model 
and show that they are the standard interactions expected from such fields. In section |9] we 
consider a more general set of defining potentials and show how Higgs masses are generated. 
We conclude with a summary and discussion. 



2 A Class of Generally Covariant Gauge Theories 

We start by giving the most compact formulation of our class of theories. This is not the 
formulation that is most suited for practical computations, but it is conceptually the simplest. 

According to our proposal, a theory that unifies gravity with gauge fields is simply the 
most general generally covariant group G gauge theory. Thus, consider a connection A 1 in 
the principal G-bundle over the spacetime manifold M. As is usual in physics literature, the 
bundle is assumed to be trivial, so the connection can be viewed as a Lie-algebra-valued one- 
form on M. The group G that we consider is a general semi-simple complex Lie group. Reality 
conditions will later need to be imposed to select a sector of the theory that corresponds to 
a particular metric signature. Note, however, that at this point there is no metric, the only 
dynamical variable of our theory is the connection A 1 . 

As we have said, the idea is to consider the most general gauge and diffeomorphism in- 
variant action that can be constructed from A 1 . The following simple construction, gener- 
alizing verbatim considerations [20] for the case of pure gravity, provides a Lagrangian with 
the required properties. Being gauge-invariant, it must only involve the curvature two- form 
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F 1 = dA 1 + (1/2) [A, A] 1 , where is the Lie-bracket and the wedge product of forms is 
assumed. Consider the 4- form F 1 A F J . This is a 4-form valued in the space of symmetric 
bilinear forms in g, the Lie-algebra of G. Choosing an arbitrary volume 4-form (vol) we can 
write F 1 AF J = (vol)Q IJ , where now Q IJ is a symmetric nxn matrix, where n = dim(g). Since 
(vol) is defined only modulo rescalings (vol) — > a(vol), so is the matrix Q IJ that under such 
rescalings transforms as Q IJ — > (l/a)Q IJ . Let us now introduce a function f(X) of symmetric 
nxn matrices X IJ with the following properties. First, the function has to be gauge-invariant: 
f(adgX) = f(X), where ad 9 is the adjoint action of the gauge group on the space of symmetric 
bilinear forms on the Lie algebra. Second, the function must be holomorphic (we work with 
complex- valued quantities). Third, and most important, the function must be homogeneous of 
degree one f(aX) = af(X). This property allows it to be applied to the 4-form F 1 A F J , with 
the result being again a 4-form. Indeed, we have f(F ! AF J ) = (vol)f(Q IJ ), and it is easy to see 
that due to the homogeneity of /(•), the resulting 4-form does not depend on which particular 
volume form (vol) is chosen. Thus, the quantity f(F I A F J ) is an invariantly-defined 4-form, 
and it can be integrated over the spacetime manifold to produce an action: 

S[A}= [ f(F J A F J ). (1) 

J M 

As we have already said, the action is complex, so later certain reality conditions will be 
imposed. 

The presented formulation (CQ) is conceptually nice, but it is very difficult to deal with in 
practice. One of the main reasons for this is that there is no natural background around which 
the theory can be expanded to produce a physically meaningful perturbation theory. This can 
be seen as follows. The first variation of the action (JTJ) is given by: 

5S = D A SA Z , (2) 

where the derivative of /(■) with respect to F 1 can be shown to make sense and is a certain 
g-valued 2-form. The second variation is given by: 

r i 8 f 8 2 f 
62s = j2dW A [6A 6A]I + dF4w DA6Al A DaSaJ > (3) 

where the second derivative of /(•) is a zero-form. Now, the most natural "vacuum" of the 
theory seems to be 

F =°- fi^ = °- W*0. (4) 

Indeed, this would indeed be a "vacuum" of the theory in the sense that the first derivative 
of the "potential" function vanishes, which then automatically satisfies the field equations 
Da(9$ / dF 1 ) = 0, and only the second derivative is non-trivial. From ([3]) we see that in this 
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case the first "mass" term is absent, and there is only the "kinetic" term for the connection. 
Thus, it seems like the perfect vacuum to expand about. However, an immediate problem with 
this vacuum is that in the absence of any background structure the second derivative in (jlj) can 
only be proportional to the Killing form g IJ , which then gives a degenerate kinetic term. So, 
there does not seem to be any way to build a meaningful perturbation theory around (jlj). 

As an aside remark we mention that the fact that the "kinetic" form in ([3]) is necessarily 
degenerate is very important for the possibility to describe gravity as a gauge theory. Indeed, 
as work [21] showed, general relativity can be put in the form (CQ) for G = SU(2) and a very 
special choice of the function /(•). At the same time, it is known to be impossible to describe 
gravity that is mediated by a spin two particle in terms of a gauge field that corresponds 
to a spin one particle. The resulution of this seeming paradox lies in the fact that the pure 
connection formulation flTJ of gravity does not allow for a well defined perturbation theory 
around Minkowski background, and so the particles that it describs are not spin one as would 
be the case in any other gauge theory. Below we shall see how the usual spin two graviton 
arises via certain "duality" trick. 

A conventional perturbative treatment for theory (CQ) is possible, but requires a rather 
strange, at least from the pure connection point of view, choice of vacuum. Thus, as we shall 
see in details in the following sections, the usual perturbative expansion around a flat metric 
corresponds in the pure connection language to an expansion around the following point: 

F 1 = 0, §j * 0. (5) 

This is a very strange point to expand the theory about, for one seems to be sitting at a point 
that is not a minimum of the "potential". However, the non-zero right-hand-side of the first 
derivative of the potential receives the interpretation of essentially the Minkowski metric, and 
a usual expansion then results. It might seem that this choice introduces a "mass" term for 
the connection, but this is not so. In fact, the second "kinetic" term is still a total derivative 
and plays no role, and there is only the "mass" term. However, as we shall see, the connection 
is no longer a natural variable in this case, and one works with a certain new two-form field B 1 
via which the connection is expressed as A 1 ~ dB 1 , so what appears as a mass term is in fact 
the usual kinetic one but for the two-form field. 

This discussion motivates introduction of a new set of dynamical fields. These are originally 
introduced via the standard "Legendre transform" trick so that integrating them out one gets 
an original action ([1]). However, one can then also integrate out the original connection field 
and obtain a theory for the new fields only. This point of view turns out to be very profitable, 
and we develop it in the next section. 
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3 Two-form field formulation 



There are at least two different ways to arrive at the new formulation. One of them is via a 
Legendre transform from (JT]), the other one by thinking about generalizations of BF theory. 

3.1 Legendre transform 

As we have already explained, we introduce a new set of fields, given by a g-valued two-form 
B 1 . The action that we would like to consider is then of BF-type and is given by: 



Here V(-) is again a G-invariant, holomorphic and homogeneous order one function of symmetric 
n x n matrices, and as such it can be applied to the 4-form B 1 A B J , with the result being again 
a 4-form. The quantity gjj is the Killing-Cartan form on q. 

Integrating out B 1 by solving its field equation: 



which is algebraic in B 1 , we get back the formulation (JTJ) with /(■) being an appropriate 
Legendre transform of V(-). However, the formulation ([6]) is much more powerful in that we 
can now choose a constant B 1 background and obtain a well-defined perturbation theory. We 
will later see how both gravity and Yang-Mills theory appear in such a perturbative expansion. 

An alternative viewpoint on the "Legendre transform" described is as follows. As we shall see 
below, the new two-form field that we have introduced is essentially the momentum canonically 
conjugate to the connection A 1 . Thus, a meaningful analogy for the relation between ([1]) and 
([6]) is the relation between Lagrangian and Hamiltonian formulation of mechanics. The former 
one uses only position variables as dynamical variables, but leads to second-derivative equations 
of motion. The later contains an independent variable - momentum, and leads to first order 
equations of motion. Thus, loosely speaking, the action (Q can be referred to as ([I]) written in 
the "Hamiltonian form" in which the momentum variable becomes an independent dynamical 
field. 

Before we proceed with an analysis of properties of the theory (J6]), we would like to present 
an alternative derivation of this action. 

3.2 Generalization of BF theory 

An alternative way to arrive at ([6]) is to consider possible ways to generalize the topological BF 
theory. For the case of G — SU(2) this was done in |22j, and here we generalize this analysis 





(6) 
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to a semi-simple Lie group. Following this reference we begin with the action 

S[A, B]=J g u B 1 A F J - U>jj B 1 A B J , (8) 

where B 1 is a two-form valued in q, F 1 is the curvature F 1 = dA 1 + \fj K A J A A K of A 1 , /J^ 
are the structure constants, and is a function (zero-form) valued in the symmetric product 
of two copies of q. At this stage this quantity is undetermined. But we should say already 
now that it is not to be thought of as an independent field to be varied with respect to, for it 
will later be fixed by Bianchi identities. Note that only the symmetric part of $ /J enters the 
action, this is why it is assumed symmetric from the beginning. Our conventions are that we 
raise and lower indices with the Killing-Cartan metric gu and its inverse g IJ . We also note 
that for a semi-simple Lie algebra we can always find a basis in which the metric is diagonal, 
i.e. gu = 8u, where Su is the Kronecker delta. 

Varying this action with respect to the connection A L and the field B L we get, respectively, 

D A B l = dB 1 + fj K A J A B K = , (9) 
F 1 = B J . (10) 

We see that the idea of the above action ansatz is to generalize BF theory in such a way that 
the equation (Q relating B and A is unchanged, while we now allow for a non-zero curvature. 
As we have already said, we do not consider a variation with respect to $ /J because we will 
later show that the Bianchi identities fix this quantity in terms of certain components of the 
two- form field B 1 . 



Let us now take the covariant exterior derivative of ffTOl) and use (|9]) together with the 
Bianchi identity D A F ! = 0. We obtain 

Da&j A B j = . (11) 

Now, the covariant exterior derivative of DaB 1 is 

D A (D A B T ) = f JK dA J A B K + f JK f* M A J A A L A B M . (12) 

Using the Jacobi identity /jjInk + fjicfki + IkiInj = 0' ^ ne equation above can be rewritten 

as 

D A (D A B I ) = fj L F J AB L , (13) 
and using equation Q and equation (TTDT) we get 

f I JL $ J K B K AB L = 0. (14) 

Let us now compute the wedge product between (jTTj) and the one- form l^B 1 , which has 
components (i^B 1 )^ = £, a B t 1 , where £ is an arbitrary vector field. We get: 

D$u Ai^ 1 AB J) = 0. (15) 
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But using it-B^ 1 A B J ^ = ^(B 1 A B J ), we can rewrite this as: 

D$ u Ai^B 1 AB J ) =0. (16) 

Let us now define the "internal" metric h IJ by means of the following relation 

B 1 A B J = h IJ (vol), (17) 
where (vol) is an arbitrary volume 4-form. We can then rewrite (Tl6l) as: 

h u D$ IJ A^(vol) = 0. (18) 
Using the definition of h IJ , we can also rewrite (CHI) as 

fj K *Lh LK = 0. (19) 

Now, computing hu DQ IJ 

hi L D& L = h IL (d& L + 2 fj K A J $ XL ) , (20) 

we can see that the second term in the right hand side vanishes because of ( 1191) and the condition 
that the Lie algebra is semi-simple. The later is used because for a semi-simple Lie algebra it 
is possible to define a Killing-Cartan metric, in our case 6u, with respect to which the object 
fuK = 8il fjK i s completely anti-symmetric. Our final result is: 

hjjd^ IJ e = 0, (21) 

which implies 

hjjd^ IJ = 0, (22) 

since £ is an arbitrary vector. 

The above equation implies that the quantities h IJ and $ /J are not independent. Let us 
define the "potential function" V := h IJ Then, 

dV = $/ j dh IJ + h u d$ IJ = $ 7J dh IJ , (23) 

where we have used (1221) . This means that: a) the potential V is only a function of h IJ , i.e., 
V = V(h IJ ) and; b) the quantities $ /J are given 

dV 

•» ' 8h» (24) 
and; c) the potential V is a homogeneous function of order one in h IJ since 

Thus, using the above definition of h IJ , and the fact that V(-) is homogeneous, we can rewrite 
the action (jSj) as 

S = j guB 1 AF J - l - V(B ! A B J ) , (26) 
which is exactly the action ([6]) we have obtained in the previous subsection. 
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3.3 Parameterizations of the potential 



As defined so far, the theory is specified by the potential function V(-). In the action it is 
applied to a 4-form, which makes things rather inconvenient in practice, since we do not have 
much experience with functions of forms. Thus, it is desirable to rewrite it as a usual function of 
a matrix. We have already discussed how to do it by introducing an auxiliary volume form, but 
it would be nice if we could avoid any arbitrariness such as that of rescalings of (vol). A possible 
way to do this is as follows. With our choice of conventions dx^ A dx v A dx p A dx a = —e^^dx 
and we have: 

B 1 A B J = ^B^B^dx" A dx u A dx p A dx a = --e^B^B^x, (27) 

where 1^ u p u is a density weight one object that does not require a metric for its definition. 
Thus, if we now define a densitized "internal metric" 

h" = \b^B j p ^ , (28) 

we can write the action as 

S[B, A] = j g u B 1 A F J + ~ V{h) d 4 x . (29) 

Thus, the potential function is now applied to an n x n matrix (densitized), and its derivatives 
can be computed via the usual partial differentiation. For example, the first variation of this 
action can be seen to be given by 



5S 



J SB 1 A ( giJ F J - 8 -^-B J \ - gjjDAB 1 A 5A J . (30) 



Indeed, the variation of the last, potential term is given by: 

2 J ah" 2 " J dh" y ' 

where the matrix of first derivatives (dV(h)/dh IJ ) is an object of density weight zero. Then, 
the field equations of our theory can be written as: 

Fi = d -m B J (32) 

dh IJ v ; 

DB 1 = dB 1 + fj K A J A B K = . (33) 

In the literature on this class of theories a different parameterization of the potential is 
sometimes used, see e.g. the original paper [17], and also the unification paper [16]. Thus, to 



11 



avoid having to take a function of forms, and/or having to work with a homogeneous function, 
one can parameterize the potential so that an ordinary function of one less variable arises. This 
can be done via a Legendre transform trick. Thus, we introduce a new variable \1/ /J that is 
required to be tracefree gu^ IJ = 0. The idea is that the matrix \1/ 7J is the tracefree part of 
the matrix of first derivatives $ 7J = (dV/dh IJ ). In other words, let us write 

$/j = */j - -gu , (34) 
n 

where is traceless. With $ /J being a function of h IJ , so is the trace part A. However, we 
can also declare A to be a function of \I/ 7J , make \1/ /J and independent variable and write the 
action in the form: 

S[B, A, *] = j giJ B<AF J -^ Ujj - ^-9u) B 1 A B J . (35) 

Varying the action with respect to \1/ /J one gets an equation for this matrix, which, after being 
solved and substituted into the action gives back (l29l) with V(-) being an appropriate Legendre 
transform of A( 1 J r ). In the formulation ( 1351) the function A( 1 J r ) is an arbitrary function of a 
tracefree matrix \I/ /J , so there is no complication of having to require V(-) to be homogeneous. 
This formulation was used in the first papers on this class of theories, but it was later realized 
that the formulation that works solely with the two- form field B 1 is more convenient. Thus, 
we do not use ( 135]) in this paper. 



4 Hamiltonian Analysis 

To exhibit the physical content of the above theory it is useful to perform the canonical analysis. 
After the 3+1 decomposition the action reads, up to an unimportant overall numerical factor: 

S = J dt J d 3 x [P^Ai - H^j , (36) 

where 

pal ._ z*c B ^ (37) 

and the Hamiltonian H is: 

- H — A[D a P aI + B^Fl - V(B ( iP a ^). (38) 

If we dealt with the pure BF theory the last "potential" term would be absent and all the 
quantities B^ a would be Lagrange multipliers. However, now the Lagrangian is not linear in 
£>o a , and, as we shall see, all but 4 of these quantities are no longer Lagrange multipliers and 
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should be solved for. The equations one obtains by varying the Lagrangian with respect to Bl a 
are: 

g«6c F / = v^p-J^ (39) 

where Vtt denotes the matrix of first partial derivatives of the function V(-) with respect to 
its arguments: 



The equations (1391) can be solved in quite a generality by finding a convenient basis in the Lie 
algebra. Thus, consider the momenta P aI . There are at least n — 3 vectors N^, a — 1, . . . , n — 3 
that are orthogonal to the momenta: 

p aI Ni = 0, Va, a. (41) 
These vectors can be chosen (uniquely up to SO(n — 3) rotations) by requiring: 

K*fl = 8*0. (42) 
We can then use the qauntities P aI , a = 1, 2, 3, N^, a = 1, . . . ,n — 3 as a basis in the Lie algebra. 
We can now decompose the quantity B^ as: 

BL = P bI B ab + NlB a a , (43) 

where B a b, B% are components of B^ in this basis. There are in total 3n components of B^ a and 
they are represented here as 9 quantities B ab as well as 3(n — 3) quantities B%. The argument 
of the function V(-) is now given by: 

B^ a P aJ) = P b{I P aJ) B ab + N { jB^P aJ] . (44) 

It is clear that this depends only on the symmetric part B ab of the components B ab . Thus, the 
anti-symmetric part of this 3x3 matrix cannot be determined from the equations (|39|) and 
thus N a in B[ a b] '■= (l/2)e a fe c A^ c remain Lagrange multipliers. It is also clear that due to the 
homogeneity of V(-) one more component of B^ a cannot be solved for. This can be chosen for 
example to be the trace part Bl a P aI , which will then play the role of the lapse function. All 
other 6 + 3(n — 3) — 1 components of Bq g can be solved for for a generic function V(-), i.e. under 
the condition that the matrix of second derivatives of V(-) is non-degenerate. We are not going 
to demonstrate this in full generality, but will verify it in the linearized theory below. 

After the quantities Bl a are solved for we substitute them into fl38l) and obtain the following 
Hamiltonian: 

-H = A^D a P al + N a P bI F I ab + NA(F, P), (45) 
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where N is the lapse function and A(F,P) is an approprite Legendre transform of V(-) that 
now becomes a function of the curvature F£ b and momentum P aI . Thus, there are n Gauss as 
well as 4 diffeomorphism constraints in the theory. It should be possible to check by an explicit 
computation that they are first class, as was done, for example for the case of G = SU(2) in 
[23J, but we shall not attempt this here, postponing such an analysis till the linearized case 
considerations. The above arguments allow a simple count of the degrees of freedom described 
by the theory: we have 3n configurational degrees of freedom minus n Gauss constraints minus 

4 diffeomoprhisms, thus leading to 2n — 4 DOF. Thus, when G = K x SU(2) the above count of 
DOF gives the right number for a gravity plus K Yang-Mills theory. For a general G one might 
suspect that the centralizer of the gravitational SU(2) describes Yang-Mills, while other part 
of the Lie algebra corresponds to some new kind of fields. Below we will unravel their nature 
by considering the linearized theory. We also note that the above count of degrees of freedom 
agrees with the one presented in [19] for the case G = SO (4). Thus, it was seen there that the 
theory describes in total 2-6 — 4 = 8 DOF, which were interpreted as those corresponding to 
2 graviton polarizations plus six new DOF. 

5 The Linearized Theory: General considerations 

As we have seen in the previous section, the mechanism that selects the gravitational SU(2) in 
G is that the momentum variable P aI provides a map from the (co-) tangent space to the spatial 
slice into g. This selects a 3-dimensional subspace in q that plays the role of the gravitational 
gauge group. Below we are going to see this mechanism at play at the level of the Lagrangian 
formulation, by studying the linearization of the action (jfJJ). In this section it will be convenient 
to introduce a certain numerical prefactor in front of this action so that the normalization of 
the graviton kinetic term in the case of gravity will come out right. Thus, we shall from now 
on consider the following action 



5.1 Kinetic term 

In this section we present some general considerations that apply to any background. We 
specialize to the Minkowski spacetime background in the next section. Let us call the first term 
in (j4"6"l) Sbf and the second "potential" term Sbb- Then, the second variation of Sbf is given 




(46) 



where i 




by: 




(47) 
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and the action linearized around Bo, A$ is obtained by evaluating this on Bq, Aq. 

As we have already mentioned, we are to view our theory as that of the two- form field B 1 , 
with the connection A 1 to be eliminated (whenever possible, see below) by solving its field 
equations. Thus, let us assume that we are given a background two-form Bq. The linearized 
connection is then to be determined from the linearized equation ([9]) that reads: 

DqSB 1 + [5A, B Y = 0, (48) 

where D is the covariant derivative with respect to the background connection Aq. Now the 
background two-form Bq is a map from the six- dimensional space of bivectors to g, and thus 
selects in g at most a 6-dimensional preferred subspace. Let us denote this subspace by 6. This 
subspace may or may not be closed under Lie brackets, but for simplicity, in this paper we shall 
assume that our background Bq is such that t is a Lie subalgebra (below we shall make an even 
stronger assumption about 6). It is then clear that the part of 6 A 1 that lies in the centralizer 
of 6 in g drops from the equation ( 1481) and cannot be solved for. As we shall later see, this will 
be the part of the group that is to describe Yang-Mills fields. The other part of 5 A 1 can in 
general be found. For this part of the connection both terms in (j4"Tj) are of the same form due 
to (j4"8"j) . and the linearized action can be written compactly as: 

S 2 S B f = 4i J 5B 1 A DoSA 1 , (49) 

where 5 A 1 has to be solved for from fj4*8l . On the other hand, for the subgroup of g that 
centralizes 6 the last term in ( 14"T1) is absent and we have: 

S 2 S BF = 8iJ 5B 1 A DoSA 1 . (50) 

Thus, our analysis of the "kinetic" term is going to be different for different parts of the Lie 
algebra. 



5.2 Potential term 



In this subsection we compute the second variation of the potential term S B b and discuss how 
it can be evaluated on a given background. We have: 

5 2 S BB = 4i / 2 ? V fy (B 5B) IJ (B 5B) KL + (5B5B) IJ , (51) 

J dh KL dh IjK J K ! dh IJ V ; V ; 

where the integration measure d 4 x is implied, and we have introduced notations 

(B 5B) IJ = \l^B^JB J p l (5B5B) IJ = \e^5B%8B J ^ (52) 
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and where the matrix of second derivatives is of density weight minus one. 

Let us now discuss how the derivatives of the potential can be computed. In general, with 
the potential function V(h) being homogeneous order one function of an n x n matrix, it can be 
reduced to a function of ratios of its invariants. A subset of invariants is obtained by considering 
traces of powers of h IJ . However, in general these are not all invariants, and other invariants 
will be introduced and discussed below in section [9j But for now, to simplify the discussion, 
let us consider a special class of potentials that only depend on the invariants obtained as the 
traces of powers of h IJ . Many aspects of our theory can be seen already for this special choice. 
Thus, consider the potential of the form: 

Trh ( Trh 2 Trh n \ 

n \(Trh) 2 (Trh) n J 

where / is now an arbitrary function of its n — 1 arguments, Tr h = gu h IJ and 

Trh? = h% a h% a h M » Ml , (54) 

for p > 2. In fact, in view of the fact that the rank of h IJ is at most six, not all the invariants 
are independent, so we could consider only 5 first arguments of /(•). Note that /(•) here is 
distinct from the function used in the action (TjQ) in the pure connection formulation of our 
theory: it is now an arbitrary function of its arguments, while this symbol in ([1]) stands for a 
homogeneous order one function. 

The parameterization given allows derivatives to be computed. Thus, the first derivative of 
the potential function with respect to h IJ is 

dh IJ n n dh IJ 

with (df/dh IJ ) given by 

df A ( d (Trh? 



dh IJ ^ p dh IJ \ (Tr h) 




Trh? 

9IJ (56) 



(Tr h)P (Tr h)P +1 

where /' is the derivative of / with respect to its argument (Tr k? /(Tr h) p ) and k?jj is 

i/j = hi Mi h Ml M2 



H.J = h IMl h% 2 h M ^j . (57) 



The second derivative of V(h) is given by: 



d 2 V(h) _ 9IJ df | qkl df | Trh d 2 f 



dh KL dh IJ n dh KL n dh IJ n dh KL dh IJ 
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with {d 2 f / dh KL dh IJ ) given by 



d 2 f V^V^f" 9 ( TrhP \ 9 ( Trh q 



dh KL dh IJ j^ 2 j^ 2 q dh IJ \{Trh)p) dh KL \{Trh) 

+ f^ 2 P dh KL dh IJ \(Trh)p ) ' 
where f" stands for the derivative of f with respect to its q argument and 

d 2 ( TrhP \ y dtfj 1 y 2 tfj 1 p 2 W££ 

dh KL d~h IJ \{Trh)P J ~(Trh)P dh KL {Tr ~h)P+ 1 9KL {Tr h)P+ 1 9 " 

p{p+l) Trh p 

+ (Trkr* 9,J9KL ' (60) 

with 

^JJ_ 1 . l M P~3 ,1.1. lM P -4 , , I U M V 

Q~ h KL 



9i(Kh L ) Ml h p 3 j+h IiK h L)Ml -- ■ -h p 4 j+ +h IMl h P \ K 9l).j ■ (61) 



With the above formulas for the fist and second derivative of the potential it is relatively easy 
to find the linearized action for any semi-simple Lie group. 



6 The G = SU(2) Case: Gravity 

As we have already mentioned, the case G = SU(2) describes (complexified) gravity theory. A 
particular choice of the potential function, see below, gives general relativity, while a general 
potential corresponds to a family of deformations of GR. In this section, as a warm-up to the 
general G case, we shall study the corresponding linearized theory. Such an analysis has already 
appeared in [23] . However, our method and goals here differ significantly from that reference. 



6.1 The metric 

To understand how G = SU(2) case can describe gravity we need to see how the spacetime 
metric described by the theory is encoded. The answer to this is very simple: there is a unique 
(conformal) metric that makes the triple B l , where i is the su{2) index, into a set of self-dual 
two-forms. This is the so-called Urbantke metric [25] 

y/=99,» ~ ^ k B; a Bi,B k p J a ^ (62) 
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that is defined modulo an overall factor. We remind the reader that at this stage all our fields 
are complex, and later reality conditions will be imposed to select physical real Lorentzian 
signature metrics. 

Alternatively, given a metric g^ u one can easily construct a "canonical" triple of self-dual 
two-forms that encode all information about g^ v . This proceeds via introducing tetrad one- 
forms 9 1 , with / = 0,1,2,3 here. One then constructs the two-forms S /J := Q 1 A 9 J and 
takes the self-dual part of S 7J with respect to IJ . The resulting two-forms are automatically 
self-dual. They can be explicitly constructed by decomposing / = (0, a) and then writing: 

E a = i9° A9 a - l -t abc d b A 9 C . (63) 

Here i = is the imaginary unit. Its presence in this formula has to do with the fact 

that self-dual quantities in a spacetime of Lorentzian signature are necessarily complex. Thus, 
even though at this stage there is no well defined signature (all quantities are complex), it is 
convenient to introduce i here so that later appropriate reality conditions are easily imposed. 
We note that "internal" Lorentz rotations of the tetrad 9 1 at the level of S a boil down to 
(complexified) SU(2) rotations of S a . 

A general sw(2)-valued two-form field B i carries more information than just that about a 
metric. Indeed, one needs 3x6 numbers to specify it, while only 10 are necessary to specify 
a metric. A very convenient description of the other components is obtained by introducing a 
metric defined by B % via (1621) and then using the "metric" self-dual two-forms (I6"3"|) as a basis 
and decomposing: 

B l = b l a Y, a . (64) 

The quantities b l a give 9 components, the metric gives 10, and the choice of "internal" frame 
for S a adds 3 more components. There is also a freedom of rescalings b\ — > Q~ 2 b l a , E a — > fi 2 E a , 
as well as freedom of SO(3) rotations acting simultaneously on S a and b l a , overall producing 18 
independent components of B l . 

When one substitutes the parameterization (1641) into the action ([6]) one finds that the fields b l a 
are non-propagating and should be integrated out. Once this is done one obtains an "effective" 
Lagrangian for the metric described by E a . Below we shall see how this works in the linearized 
theory. However, we first need to choose a background. 

6.2 Minkowski background 

The Minkowski background is described in our framework by a collection of metric two-forms 
(1631) constructed from the Minkowski metric. Thus, we choose an arbitrary time plus space 
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split and write: 



Eg = idt A dx a - \e abc dx b A dx c , (65) 



where dt, dx a , a = 1, 2, 3 form a tetrad for the Minkowski metric ds 2 = —dt 2 + 'Yl, a {dx a ) 2 . Our 
two-form field background is then chosen to be 

B l = <5*Eg, (66) 

where 5* is an arbitrary SO (3) matrix that for simplicity can be chosen to be the identity 
matrix. 

In what follows we will also need a triple of anti-self dual metric forms that, together with 
(1631 form a basis in the space of two-forms. A convenient choice is given by: 

So = idt A dx a + ^e abc dx b A dx c . (67) 

The following formulas, which can be shown to follow directly from definitions (165]) and 
( 1671) . are going to be very useful 

(68) 
(69) 
(70) 
(71) 
(72) 

where r]^ is the Minkowski metric. We are going to refer to them as the algebra of E's. 

The first of the relations above, namely (168]) . is central, for all others (apart from (172]) ) can be 
derived from it. It is useful to develop some basis-independent understanding of this relation. 
We are working with the Lie algebra su(2) and are considering a basis X a in it in which 
the structure constants read [X a ,X 6 ] = e abc X c . This is the basis given by X a = — (i/2)<j a , 
where a a are Pauli matrices. The metric g ab = 5 ab on the Lie algebra can be obtained as 
g ab = — 2TT(X a X b ). Then (168]) can be understood as follows: the product of two E's is given 
by minus the metric plus the structure constants times E. We will see that in this form the 
relations (168]) persist to any basis in su(2). 



6.3 Linearized action 

We are now going to linearize the G = SU(2) theory around the background (166]) . Thus, we 
take: 

B i = B l + b i . (73) 





= -£ a V + 6 abc E^, 




= 45 ab } 


abc y^a spba spc 
e ^OAtcr^O A^OA/i 


= -4!, 


abc y-ia y-ifr spdva 


= -25 cd Vlip , 


y^a spa 


^Ifip^lua ^]fi<T^]up ^p.vpa i 
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As we have already discussed, to linearize the kinetic BF term of the action we need to solve 
for the linearized connection if we can. This is certainly possible for the case at hand, as we 
shall now see. 

If we denote the linearized connection by a 1 we have to solve the following system of equa- 
tions 

dtf + e) k a j A £>q = , (74) 

where we have used the fact that the background connection is zero. It is convenient at this 
stage to replace all i-indices by a-ones, which we can do using the background object 5 l a that 
provides such an identification. We can now use the self-duality e^^Y^^ = liY?^ of the 
background to rewrite this equation as 

1 



-e^d u b a pa + e abc a h u Y^ u = 0. (75) 

We now multiply this equation by EL, and use the identity (ITT]) to get: 

= a K a ^ up(7 d u b b pa , or a% = ±Z b 0/3 "Eg^^", (76) 

where we have introduced a compact notation: 

(db b y := e^ pa d u b b pa (77) 
for a multiple of the Hodge dual of the exterior derivative of the perturbation two-form. 



The BF part of the linearized action was obtained in ff49|) . We need to divide the second 
variation given in this formula by 2 to get the correct action quadratic in the perturbation. 
Thus, we have: 



S® F = 2i J b a A da a = -ij a a ^db a Y, (78) 

where we have written everything in index notations and integrated by parts to put the deriva- 
tive on b^v, and used the definition (1771) . Now substituting (ITS"]) we get: 

s [ b ] f = \J r"K a ,{db b Y^ Pv {db a y. (79) 

Let us now linearize the potential term. For this we need to know the background h l i as well 
as the matrices of first and second derivatives for the background. Using (1651) is easy to see that 
h l Q = 2i5 l: ' . Since the background volume form is just the identity we can now safely remove 
the density weight symbol from the matrix h^. Also, as before, let us replace all z-indices by 
a-indices using 5*. Using (155]) and the fact that the first derivatives (df/dh ab ) vanish on this 
background we immediately get: 



dV 



Q h ab 



= ^ fo, (80) 

h 6 
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where fo is the background value of the function / in the parameterization (1531) . It is not 
hard to see that this value plays the role of the cosmological constant of the theory, so in our 
Minkowski background it is necessarily zero by the background field equations. The matrix of 
second derivatives of the potential is easily evaluated using (1551) and we find: 



d 2 V 



dh cd dh ab 

where we have introduced: 



ho 



g;=E W-i) , (g2) 

p=2,3 

This is a constant of dimensions of the cosmological constant 1/L 2 . It is going to play a role of 
a parameter determining the strength of gravity modifications. 



We can now write the linearized potential term (loTl) . We must divide it by two to get the 
correct action for the perturbation. This gives: 

S { b ] b = "f / (6a(c6 d)b - i^c) (Er&W W%) ■ (83) 

Note that the tensor in brackets here is just the projector on the tracefree part. This fact will 
be important in our Hamiltonian analysis below. Our total linearized action is thus (1791) plus 



6.4 Symmetries 

The quadratic form obtained above is degenerate, and its degenerate directions correspond 
to the symmetries of the theory. These are not hard to write down. An obvious symmetry 
is that under (complexified) SO (3) rotations of the fields. Considering an infinitesimal gauge 
transformation of the background T/q we find that the action must be invariant under the 
following set of transformations: 

<U*„ = e abc co%^ (84) 

where u a are infinitesimal generators of the transformation. It is clear that (183)) is invariant 
since it involves only the aft-symmetric part of (Sq^o^), and the transformation (1541) affects 
the anti-symmetric part. Let us check the invariance of the kinetic term ( 1791) . We have the 
following expression for the variation: 

\ J V a ^ a 0a ^b b yZ b opu (dbr- (85) 
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Substituting here the expression for the variation we find: 

^S^WyX^ = 2i?f ^ a / CC V c X^XU = m»\ (86) 

where we have used the self-duality of £quk and applied the identity ( JTT1) once. Substituting 
this to (1851) and integrating by parts to move the derivative from u a to b a we get under the 
integral & vpo d^djf = 0, since the partial derivatives commute. This proves the invariance 
under gauge transformations. 

Another set of symmetries of the action is that of diffeomorphisms. These are given by: 

5(b a = di^ a , (87) 

where l% is the operator of interior product with a vector field It is not hard to compute 
this explicitly in terms of derivatives of the components of the vector field. However, we do not 
need all the details of this two-form. Indeed, let us first note that the first "kinetic" term of 
the action is in fact invariant under a larger symmetry: 

S v b a = d V a , (88) 

where r] a is an arbitrary Lie-algebra valued one-form. Indeed, this is obvious given that the 
kinetic term is constructed from the components of the 3-form db a given by the exterior deriva- 
tive of the perturbation two- form. Thus, ( l88l) indeed leaves the kinetic term invariant. Then, 
since ( 1H7I) is of the form ( IHBl with i] a = we have the invariance of the first term. To see 
that the potential term (|83|) is invariant we should simply show that the symmetric tracefree 
part of the matrix (So^&) ab is zero. Let us compute the symmetric part explicitly. We have: 

E^d^tu = rt p£ P , (89) 

where we have used (|68|) . Thus, there is only the trace symmetric part, so the part that enters 
into the variation of the action (1851) is zero. This proves the invariance under diffeomorphisms. 
Note that the second "potential" term is not invariant under all transformations (IHBl . since for 
such a transformation that is not a diffeomorphism the matrix (T l o5 ri b) ab contains a non-trivial 
symmetric tracefree part, as can be explicitly checked. 

We will see that these are the only symmetries when we perform the Hamiltonian analysis. 
However, before we do this, let us show how the usual linearized GR appears from our theory. 

6.5 Relation to GR 

In this subsection we would like to describe how general relativity (linearized) with its usual 
gravitons appears from the linearized Lagrangian described above. We shall see that to get 
GR we must take the limit when the "mass" parameter g for the components (Sq6)^, where 
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tf stands for the tracefree part, is sent to infinity. Indeed, the potential part (15311 depends 
precisely on these components, and when the parameter g is sent to infinity these components 
are effectively set to zero. We shall now see that this gives GR. 

It is not hard to show that in general the tracefree part := h^ u — {1/ A) , q^ v h p p of the 
metric perturbation defined via = r/^ + corresponds in our language of two forms 
to the anti-self-dual part of the two-form perturbation: 

(ay- = s o[/^ P - (so) 

The fact that this two-form is anti-self- dual can be easily checked by contracting it with Eq MW 
and using the algebra ( J68l) . The result is zero, as appropriate for an anti-self-dual tow-form. 
In addition to (I9~0l) there is in general also the self-dual part of the two-form perturbation. 
However, in the limit g — > oo all but the trace part of this gets set to zero by the potential 
term. The trace part, on the other hand, is proportional to the trace part rfh^v of the metric 
perturbation. To simplify the analysis it is convenient to set this to zero rf^h^ = 0. This is 
allowed since in pure gravity the trace of the perturbation does not propagate. Then (J90l) is 
the complete two-form perturbation, and we can drop the tf symbol. 

To simplify the analysis further, instead of deriving the full linearized action for the metric 
perturbation h^, let us work in the gauge where the perturbation is transverse d^h^ = 0. Let 
us then compute the quantity (<% a ) M in this gauge. Using anti-self-duality of b^ v given by (1901) 
we have: 

e» vp °d v b a pa = -2id v b atlv . (91) 
Substituting here the explicit expression (190]) and using the transverse gauge condition we get: 

{db a Y = YL a Q lJp d u h p l . (92) 
We can now substitute this into the action ( 1791 to get: 

S (2) =-\J V a "KaA PCTd p h ^UuK l5 d,K (93) 

= -\j n aP {5X - * 8 Jl - K^W^K - Sffi - ie^HW*. 

where we have used (1721 to get the second line. We can now contract the indices and take into 
account the tracefree as well as the transverse condition on h^. We get the following simple 
action as the result: 

S {2) = ~f d„h pa dW, (94) 

which is the correctly normalized transverse traceless graviton action. Note that in the passage 
to GR we have secretly assumed that in (1901) is a real metric perturbation. Below we will 
see how to impose the reality conditions on our theory that this comes out. Also note that the 
sign in front of ( 1941) is correct for our choice of the signature being (—,+,+, +). 
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6.6 Hamiltonian analysis of the linearized theory 

For a finite g our theory describes a deformation of GR. Since not all components of the two- 
form perturbation 6" are dynamical, the nature of this deformation is most clearly seen in the 
Hamiltonian framework. This is what this subsection is about. 

We note that the outcome of this rather technical subsection is that at "low" energies 
E 2 <C g the modification can be ignored and one can safely work with the usual linearized GR. 
Thus, it may be advisable to skip this subsection on the first reading. Let us start by analyzing 
the kinetic BF-part. 



Kinetic term. Expanding the product of two E-matrices in (!79l) using fl68l) we can write the 
linearized Lagrangian density for the BF-part as 

C BF = ±(db a r(db b y (e abc i: c 0flI/ + 6 ab V , u ) . (95) 

Let us now perform the space plus time decomposition. Thus, we split the spacetime index as 
fi = (0,a), where a = 1,2,3. Note that we have denoted the spatial index by the same lower 
case Latin letter from the beginning of the alphabet that we are already using to denote the 
internal su(2) index. This is allowed since we can use spatial projection of the £q mj , two-form 
to provide such an identification. Thus, from fl63l) we have: 

K bc = -t\a (96) 

and 

Kob = < ■ (97) 



Let us now use these simple relations to obtain the space plus time decomposition of the 
Lagrangian. First, we need to know components of the (db a y vector. The time component is 
given by: 

(db a f = e obcd d b b a cd = -d b t ab , (98) 
where our conventions are e 0abc = —e abc and we have introduced: 

The spatial component of (<% a ) M is given by: 

(db a ) b = e Wcd d b a cd + 2e bc0d d c b a 0d = d t ab - 2e bcd d c b a 0d . (100) 



Now, the Lagrangian (19"oT) is given by: 

Cbf = -\{db a f{db a f + l(db a )°(db b ) d e abc i; c Qd + \(db a ) e (db b ) f (e abc E c ef + S ab 5 ef ). (101) 
4 2 4 1 
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Substituting the above expressions we get: 



Cbf = - \d b t ab d c t ac - l -d d t ad {d a t bc - 2e^d e b b 0f )e abc (102) 
-hd t ae - 2e emn d m b a 0n )(d t b f - 2e^d p b b q )(e abc e c ef - 5 a % f ). 



Our fields are now therefore b^ b and t ab . There will also be another, potential part to 
this Lagrangian, but it does not contain time derivatives, so the conjugate momenta can be 
determined already at this stage. Thus, it is clear that the field b$ b is non-dynamical since the 
Lagrangian does not depend on its time derivatives. The momentum conjugate to t ab , on the 
other hand, is given by: 



7T 



ab ._ d ^BF _ _?_abc f) ,cd _ hx+ef 



- -^ ahc d d f d - -(d t £ f - 2e^d p b e 0q )(e^ - S"&). (103) 



It is not hard to check that the momentum variable is simply related to the spatial projection 
of the connection ( 1761) as: 

K = ~2< (104) 



To rewrite the Lagrangian in the Hamiltonian form one must solve for the velocities dot ab 
in terms of the momenta 7r afe . However, it is clear that not all the velocities can be solved for 
- there are constraints. A subset of these constraints is given by the \x = component of the 
( 1751) equation that, when written in terms of 7r a6 , becomes: 

ga ._ e abc n bc + ^ab = Q ( 1Q5 ) 

These are primary constraints that must be added to the Hamiltonian with Lagrange multipli- 
ers. 

Thus, the expression for velocities in terms of momenta will contain undetermined functions. 
These functions are simply the aJJ components of the connection, as well as (at this stage 
undetermined) 6q 6 components of the two-form field. The expression for velocities is given by 
the spatial components of equation (J75|) . After some algebra it gives: 

d t ab = 2e bef d e b a 0f - 2e abc a c - e aed e dbf 7i ef . (106) 



Let us now obtain a slightly more convenient expression for the Lagrangian. Indeed, recall 
that using the compatibility equation between the connection and the two-form perturbation, 
we could have chosen to write our linearized action (1751) as 



SP F = -2i J e abc ^ A a b A a c = -2 J E a " L 't abc a\a c v . (107) 
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Introducing the time plus space split and writing the result in terms of the momentum variable 
(I104p we get the following Lagrangian: 

r c\abc„ab„c ^ net nbcbecf f T no\ 

Lbf — ~2t 7r clq — — e e n 7r j . (lUo J 

We can now easily find the BF-part of the Hamiltonian: 

Ubf = 7r ab d t ab - C BF = 2n ab e bef d e b a 0f - ^e aei 'e abc n be n cj ". (109) 

We need to add to this the primary constraints (11051) with Lagrange multipliers. Thus, the 
total Hamiltonian coming from the BF part of the action is 

H totai = 27r ab e bef d e b% f - ~e aef e abc ir be Tr cf + u a Q a . (110) 
1 2 

This is, of course, the standard result for the linearized BF Hamiltonian. If not for the potential 
term, the Hamiltonian would be a sum of terms generating the topological constraint c^vr^ = 
and the Gauss constraint (11051) . Let us now consider the other BB part of the Lagrangian. 

Potential part. We can rewrite the linearized Lagrangian density for the BB part fl83l) as 

= (^) tf (* s 5 w ) tf . (in) 

where tf stands for the tracefree parts of the matrices. Splitting the space and time indices 
gives: 

(bil^ u ) t{ =-(2ib^ b) +t^) t , (112) 



and so 



jC bb = -- ( 2ib {ab) + t( ab A ( 2ib ( f + . (113) 

2 V /tf \ /tf 



Analysis of the constraints. Thus, the total linearized Hamiltonian density H = W'^f 1 — Lbb 

is given by 

H = 2li ab e bef d e b a f - ^fe^be^cf + ^aga + 9 r 2ib (ab) + f ^(ab) + j (a6) \ _ 

2 2 V /tf \ / tf 

It is now clear that only the anti-symmetric part and trace parts of 6 a g remain Lagrange 
multipliers in the full theory. These are the generators of the diffeomorphisms. The other part 
of b a Q, namely the symmetric traceless is clearly non-dynamical and should be solved for from 
its field equations. Varying the Hamiltonian with respect to this symmetric tracefree part we 
get 

(2ib iab) +t^ b A = i -(e e ^d e 7r b A . (114) 
V /tf g \ 1 / tf 
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Now writing: 

b ab Q = iN5 ab + ^e abc N c + (b ( f) tf (115) 
and substituting the symmetric tracefree part from (11141) we get the following Hamiltonian 

H = -2Nie abc d a 7t bc - 2d [a ^ } N b + u a G a (116) 

_ laef^bc^cf + j f e ef(a de7r b)\ ^ a b)^ + }_ Lef{aQj>)\ ^(ag^ _ 

The reason why we introduced a factor of i in front of the lapse function will become clear 
below. One can recognize in the first line the usual Hamiltonian, diffeomorphism and Gauss 
linearized constraints of Ashtekar's Hamiltonian formulation of general relativity [15]. The first 
two terms in the second line comprise the Hamiltonian. Finally, the last term is due to the 
modification and goes away in the limit g — ► oo. 

It is not hard to show that the reduced phase space for the above system is obtained by 
considering 7r afc , t ab that are symmetric, traceless and transverse d a 7r ab = 0, d a t ab = 0. On such 
configurations the matrix e e * a d e irf b is automatically symmetric traceless and transverse. The 
reduced phase space Hamiltonian density is then given by: 

H phys = _( 7r a6)2 + ^fage^ab + ,ga^2 
2 2g 

where we have integrated by parts and put the derivative on t ab in the second term. This 
Hamiltonian is complex, so we need to discuss the reality conditions. 

Reality conditions. So far our discussion was in terms of complex-valued fields. Thus, the 
reduced phase space obtained above after imposing the constraints and quotienting by their 
action was complex dimension 2 + 2. Reality conditions need to be imposed to select the 
physical phase space corresponding to Lorentzian signature gravity. 

In the case of GR that corresponds to g — > oo the reality condition could be guessed from 
the form of the Hamiltonian H 1 1 7f) . Indeed, we can write it as: 

U P g V r s = - (7c ab + ie efa d e t fb ) 2 + ^(d a t bc ) 2 . (118) 

Thus, it is clear that we just need to require t ab and ir ab + ie e ^ a d e t^ b to be real. This procedure, 
however, does not work for the full Hamiltonian because of the last term in (11171) . 

Let us now note that the last term in (I117p . when written in momentum space behaves as 
E 2 /M 2 , where E is the energy and M 2 = g is the modification parameter. Thus, for energies 
E M the modification term is much smaller than the term it 2 and can be dropped. It is 
natural to expect that gravity is only modified close to the Planck scale, so it is natural to 
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expect M 2 « M 2 , where M p is the Planck mass. With this assumption the last term in ( 11 17ft is 
unimportant for "ordinary" energies and can be dropped. Thus, if we are to work at energies 
much smaller than the Planck scales ones then we do not need to go beyond GR described by 
the first two terms in (j!17p . 

The above discussion shows that a discussion of the reality conditions for the full Hamil- 
tonian (11171) . even though possible and necessary if one is interested in the behavior of the 
theory close to the Planck scale, is not needed if one only wants to work for with much smaller 
energies. For this reason, and in order not to distract the reader from the main line of the 
argument, a somewhat technical reality conditions discussion for the full theory is placed in the 
Appendix. 

Now that we understood how the simplest case G = SU(2) gives rise to gravity we can apply 
the same procedure to more interesting cases of a larger gauge group. We consider the example 
of SU(3) that well illustrates the general pattern. 



7 The G = SU(3) Case: Gravity- Maxwell system 

In this section we perform an analysis analogous to that in the previous section but taking a 
larger gauge group. As before, we first consider the complex theory, and only at the end impose 
the reality conditions. Let us start by reviewing some basic facts about the su(3) Lie algebra. 



7.1 Lie algebra of SU(3) 

The standard matrix representation of the Lie algebra of SU(3) consist of all traceless anti- 
hermitian 3x3 complex matrices. The standard basis for su(3) space is given by the imaginary 
unit times a generalization of Pauli matrices, known as Gell-Mann matrices. These hermitian 
matrices are given by: 



Ai = 1 , A 2 = i00 , A 



A 4 = , As = , A 






A 7 = -i , A 8 = -= 1 . (119) 
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Table 1: Commutators between T + , T_,T Z , V+, VI, W+, W_, F. 



However, in our computations the Cartan-Weyl basis is going to be more convenient. Let 
us recall that in the Cartan-Weyl formalism one starts with the maximally commuting Cartan 
subalgebra, which in our case is spanned by two elements A 3 , A 8 . One then selects basis vectors 
that are eigenstates of the elements of Cartan under the adjoint action. This leads to the 
following basis, see [26].[27] 



-L(T x ±iT y ) V ± =-L(V x ±iV y ) w ±=^ 



where T x = |Ai, T y = ~A2, V x = ~A4, V y = ^X^, W x = ^Xq and W y = 7^7. Then the Cartan 
subalgebra is Hi = Span(T z , Y), and the commutator between any of the H^s and the rest 
of the elements of the basis E a , E a = {T + , T_, T z , V + , VI, W + , W_}, is a multiple of E a , i.e. 
[Hi,E a ] = a.iE a . One considers the a^s, for i = 1,2, as the components of a vector, called 
a root of the system. In this case we have six roots, i.e. {1,0}, { — 1,0}, {-, ^}, {— |, — 
{— |, 2 }> {|> — 2 }■ ^he brackets between elements of this basis are given in Table [TJ We 
also need to know the metric gjj = — 2Tr(T/Tj) in this basis. It is given in Table [2j 



7.2 Background 

Let us now discuss how a background to expand around can be chosen. A background two- 
form field Bq is a map from the space of bivectors, which is 6-dimensional, to the Lie algebra in 
question. Thus, its image is at most 6-dimensional subspace in su(3). There are many different 
subspaces one can consider. In this paper we study the simplest possibility. Thus, we choose 
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Table 2: Components for the internal metric in the base {T + , T_,T Z , V+, VI, W+, W-, Y}. 



Bq such that the image of the space of 2-forms that it produces in su(3) is 3-dimensional. 
Moreover, we choose this image to be an su(2) Lie sub-algebra. Even further, we choose this 
sub-algebra to be that spanned by {T + , T_, T z }. Clearly, this is not the only su{2) sub-algebra 
in su{3). Other possibilities include {V + , VI, \ (y/3Y + T z )} and {W+, \ (y^Y - T z )}. In 
this paper we do not study these different possibilities, leaving a more thorough investigation 
to further research. We believe that the example we choose to study is sufficiently illustrating. 

Thus, our background is essentially the same as the one we considered in the previous 
section. This is motivated by our desire to have the usual gravity theory arising as the part 
of the larger theory we are now considering. Since in the general Lie algebra context it is 
convenient to work with the Cartan-Weyl basis, we need to change the basis of basic two-forms 
(I65p as well. This can be worked out as follows. In the previous section we were using a 
basis in the Lie algebra in which the structure constants were given by e a b c - If we denote the 
corresponding generators by X a then [X a ,X fe ] = e ahc X c . On the other hand, for generators T a 
used in (11201) we have [T a ,Tj,] = ie a b c T c . The relation between these two bases is X a = —\T a . 
We can then define a new set of self-dual two-forms E*, T, z via: 

e = sax " = S+T + + S ~ T - + SZ ^- ( 121 ) 

a=l,2,3 

This gives 

E+ = ^(E 1 -iE 2 ) E- = ^ (E 1 + iE 2 ) E^ = -iE 3 . (122) 

y/2 V ' V2 K ' 

The sw(3)-valued two-form E is our background to expand about. 
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7.3 Linearization: Kinetic term 



As before, the first step of the linearization procedure is to solve for those components of the 
connection for which this is possible. As we have discussed in section [5j this is in general 
possible for the components of the connection in the directions in the Lie algebra that do not 
commute with the directions spanned by the background two-forms. In our case these are the 
directions spanned by T±,T Z and V±,W±. We already know how to solve for the connection 
components in the directions T±,T Z . Indeed, the solution is given by ( !76l) which we just have 
to rewrite in the different basis. It is, however, more practical to solve the equations once more 
by working in the different basis from the very beginning. 

The su{2) part. The su(2) sector equations in the Cartan-Weyl basis are: 

db + + a z A £ + - a + A £ z = 0, 

db~ + a~ A 57 - a z A £~ = 0, (123) 
db z + a + A £~ - a' A E+ = 0. 

We rewrite them in spacetime notations, take the Hodge dual, and use the self-duality of the 
S ± , S z matrices to get: 

^{db + Y + <£ +/u/ - a+E*"" = 0, 

h.db-Y + a-H z ^ -alTT^ = 0, (124) 
^{db z Y + a+E-* 4 " - a;S + ^ = 0, 

where the notation is, as before (<%) M = e^ vpa d v b pa . We now need the algebra of the new 
E-matrices. It can be worked out from the relations f)122p and the algebra (|68|) . We get: 

£^£-=r^, E+E+^ = 0, E^£-- = 0. (125) 

For purposes of the calculation it is very convenient to rewrite these relations in the schematic 
form, by viewing them as matrix algebra. Our matrix multiplication convention for the two- 
forms is (XY)^ = XPY p v . We have: 

£+£- = ?? + £*, £*£+ = £+, £ 2 £~ = -E~ , 

£ 2 £ 2 = r/, E + E + =0, E"E" = 0. (126) 

This is precisely the relations fl68l) . just written in terms of metric and the structure constants 
on su(2) for a different basis. 
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In matrix product conventions, the equations (11241) take the following transparent form: 

S z a+ = 0, 

Z-a z = 0, (127) 
£ + a~ = 0, 

where the convention is that the second spacetime index of £ is contracted with the spacetime 
index of a. 

We can now solve (11271) by using the algebra (I126p . To this end we multiply the first 
equation by S + and the second one by £~. This leads to two equations involving only a ± but 
not a z . We can obtain another two equations of the same sort by multiplying the last equation 
in (11241) by H ± . Then adding-subtracting the resulting equations we get: 

a+ = ~ (£-£+(96+) + £+(<%*)) , a- = ~ (£+£-(0&-) - E~{db z )) . (128) 

To obtain the last component of the connection we multiply the first equation in (I127p by 
£~ and second by X + , and then subtract the resulting equations. We find E~a + — £ + a~ = 
-(l/2i)(<% z ) using (1125)) . We get: 

a z = -i {{db z ) + £-(<%+) - E + {db~)) . (129) 

It is now easy to write the su(2) part of the linearized BF part of the action. Using the 
metric components given in Table [2} from ( l78l) we have: 

5^ 2) = ~J (db + ) (E + E-(db-) - £-(<%*)) + (db~) (S-S+(96+) + S + (96 2 )) (130) 

+{db z ) ((db z ) + £"(<%+) - E+{db-)) , 
where again our convenient schematic form of the notation is used. This is simplified to give: 

= ~J (db + )(ri + Z z )(db-) + (<%-)£+(<% z ) - {db + )E-{db z ) + ^(db z )(db z ). (131) 

We could now use this as the starting point of the Hamiltonian analysis similar to the one in 
the previous section. However, it is clear that its results are basis-independent, so we do not 
need to repeat it. Still, the above considerations are quite useful as a warm-up for the more 
involved analysis that now follows. 

The part that does not commute with su(2). Let us denote the four directions V±,W± 
collectively by index a = 4,5, 6, 7. We have to solve the following system of equations: 

db a + /1/Ar = 0, (132) 



h,db + ) + Y, + a z - 
L(db-)+Z*a-- 
hdb z ) + Y,~a + ~ 
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where the terms f a ab a a A E 6 are absent since the corresponding structure constants are zero. 
Explicitly, using tabled] we have: 



db 4 - 




--a 4 AE z 
2 


= 0, 


(133) 


db 5 + 


1 7 

a A E~ 

V2 


+ -a 5 AE 2 
2 


= 0, 


(134) 


db 6 - 


\/2 


+ i a 6 A E* 


= o, 


(135) 


db 7 + 


v/2 




= 0. 


(136) 



We can solve this system using the same technology that we used above for the su(2) sector. 
Thus, we take the Hodge dual of the above equations, use the self-duality of the E's, and rewrite 
everything in the schematic matrix form. We get: 

-(db 4 ) - 4= S + a 6 - - E 2 a 4 = , 
2i v ; v/2 2 

^ 6 )-^S-a 4 + ls¥ = 0, (137) 

-(96 7 ) + ^E + a 5 --EV = 0. 
2i V ; v/2 2 

We can now manipulate these equations using the algebra (11261) . Thus, let us multiply the 
third equation by y/2T. + and subtract the result from the first equation. This gives: 

i(96 4 ) - ^S+(«% 6 ) + (v + = 0. (138) 

It is now easy to find a 4 by noting that (rj + (1/2)E 2 ) -1 = (4/3) (77 — (1/2)E 2 ). Thus, we have: 

a 4 = 1 (^v / 2E + (96 6 ) - (2 V - S*)(«% 4 )) . (139) 

Similarly, we multiply the last equation by v^E" and add it to the second equation. Mul- 
tiplying then by the inverse of (77 — (1/2)E 2 ) we get: 



a 5 



~ (V2Z-(db 7 ) + (277 + S 2 )(<% 5 )) . (140) 



To find a 6 we multiply the first equation by y/2T, and subtract the result from the third 
equation. We then multiply the result by the inverse of (77 — (1/2)E 2 ). We get: 

a 6 = - ( V2Z-(db 4 ) - (2rj + Z z )(db 6 )) . (141) 

61 \ J 
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Finally, to find a 7 we multiply the second equation by \/2Y> + and add the result to the last 
equation. Multiplying the result by the inverse of (77 + (1/2)T> Z ) we get: 

a 7 = -1 (^v / 2S+(96 5 ) + (277 - Z z ){db 7 )) ■ (142) 

We should now substitute the above results into the relevant part of the action. This is 
again obtained from (1781) by taking into account the expression for the metric. We shall refer 
to this part of the action as "Higgs" in view of its interpretation to be developed later. We 
have: 

Sf F " s = i J a 4 (8b 5 ) + a 5 (<% 4 ) + a\db 7 ) + a 7 (<% 6 ), (143) 

where we took into account and extra minus sign that comes from the metric. Substituting 
here the above connections, we get, after some simple algebra: 

S™ 998 =H J V2(db 5 )E + (db e ) - V2{db 4 )^(db 7 ) (144) 
-(db 4 )(2r] + £ z )(<% 5 ) - (<% 6 )(2r? - Y, z ){db 7 ). 
A more illuminating way to write this action is by introducing two two-component fields: 

s) (:>)■ 

It is not hard to see that this split of the "Higgs" sector part of the Lie algebra is just the split 
into two irreducible representation spaces with respect to the action of the gravitational SU(2). 
In terms of these columns the above action takes the following form: 

Below we will use this action as the starting point for an analysis that will eventually exhibit 
the physical DOF propagating in this sector. 

Centralizer U(l) part. We cannot solve for the components of the connection in the part 
that commutes with su(2). In our case this is the direction Y of the Lie algebra. We shall refer 
to this part of the action as "YM". Thus, the action remains of BF type: 

S™ = _4i J b 8 Ada 8 , (147) 

where the extra minus sign is the one in the metric. 
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7.4 Linearization: Potential term 



As in the SU(2) case our background internal metric h 1 ^ is just 2\g ah in the su(2) directions 
and zero in all other directions. Since the background metric is flat we shall drop the tilde from 
h IJ in this section. We compute the matrix of first derivatives of the potential using (1551) . We 
get: 



dV 



Q h ab 

dV 



fo 

Z ~a~ 9abi 



dh aa 
dV 



dh a P 



fo 



V 



p=2 



3 p- 



9a(3 



(148) 
(149) 

(150) 



Here fo, (fp)o are the value of the function and its derivatives at the background, and index 
a stands for all directions in the Lie algebra that are not in su(2). The quantity fo can be 
identified with a multiple of the cosmo logical constant. More specifically: 



a = -3*. 



Let us also define another constant of dimensions 1/L 2 : 



p_ 

3P~ 



Then we have: 



dV 



dh a(3 



•(A/3 + n)g a p 



(151) 



(152) 



(153) 



The sum here and in the previous formula is taken over p — 2, ... ,6, because the function / can 
at most depend on 5 ratios of 6 invariants of the matrix h IJ . It has at most only 6 independent 
invariants since it is constructed from the map B 1 that has the rank at most six. Since we 
want to work with the Minkowski spacetime background we should set A = 0, which we do in 
what follows. 

We now need to compute the matrix of second derivatives. Let us first obtain its su(2) part. 
Using fl58l) we get: 

d 2 V 



9_ 

dh cd dh ab 2i \ Sa ^ d ^ b ^9ab9cd 



where we have defined: 



1 6 



pip - 1) 



p=2 



(154) 



(155) 
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As in the SU(2) case this constant is going to measure strength of gravity modifications. Both 
n and g constants have dimensions of 1/L 2 and are, in general, independent parameters of our 
linearized theory, related to first derivatives {fp)o of the function / of the ratios. 

Let us now compute the matrix of second derivatives in its part not in su(2). We only need 
its mixed components aa and b/3. The computation is easy and using (|58l) we get: 

d 2 V 

= T7 QabQaP- (156) 
o 41 



dh aa dh b P 



We note that in this computation only one of the terms in (1611) survives, and this is the reason 
why it is the constant k that appears in this formula. 

We can now compute all the potential parts. We use ( 15T1) which we have to divide by two 
to get the correct quadratic action. For the su{2) gravitational part the result is unchanged 
from that in the previous section and we have: 

S 9 B r B =~J (ga(cg d)b - \g ab g c )j (sr&y (£o P X) • ( 157 ) 

The "Higgs" and "YM" parts of the potential term are both given by: 

S m 99 s-YM = _kJ g ah g a ^b^^%%) + 2i^e^&^, (158) 

so the indices a, (3 here take values 4, 5, 6, 7, 8. We can further simplify this using fl72|) . We get: 

S m 9gs -YM = _ k J gaph ^ P °p- vpa . (159) 



where 



P = \ \ + ^ ( 16 °) 

is the anti-self-dual projector. 
7.5 Symmetries 

We have seen that the su(2) sector of the theory is completely unchanged from what we have 
obtained in the G = SU(2) case. One can moreover see that diffeomorphisms still act only 
within this sector. Indeed, the action of a diffeomorphism in the direction of a vector field 
£ M is still given by ( 1H71) and only changes the su(2) part of the two-form field. Similarly, the 
SU(2) gauge transformations act only on the su(2) sector. Thus, the gravity story that we have 
considered in the previous section is unchanged. 
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Let us now consider what happens in directions not in su(2). Let us first consider the "Higgs" 
sector spanned by V±, W±. A gauge transformation with the gauge parameter to valued in this 
sector acts as d^b = [u, £]. In components this reads: 



The remaining part of the Lie algebra is that spanned by Y . The corresponding gauge 
transformation has no effect on the two-form field b 8 (nor on b a , a = 4,5, 6, 7) since it commutes 
with the background. However this gauge transformation does act on the connection a 8 by the 
usual U(l) gauge transformation a 8 — > a 8 +du 8 . The kinetic part (I147P clearly remains invariant, 
and the potential part is invariant since it only depends on b 8 that does not transform. 

7.6 Low energy limit of the "Higgs" sector 

Our analysis of the "YM" sector presented below will show that the parameter k that appeared 
in the " Higgs- YM" part of the potential (11581) must be taken to be of the order M£, where M p 
is the Planck mass. This will be follow from the fact that the YM coupling constant should 
be of order unity in a realistic unification scheme, which then immediately implies k ~ M p . 
Another way to reach the same conclusion is to note that M p is the only scale in our problem, 
so all dimensionful quantities must be of the Planck size, see more on this in the last discussion 
section. If this is the case then the role of the potential term (11581) for the "Higgs" sector 
is to make the anti-self- dual components of the two-forms 6" "infinitely massive" and thus 
effectively set them to zero. This is quite similar to what happened in the gravitational sector 
in the limit g — > oo with the 6"j? components. Thus, we see that in the low energy limit £ 2 < k 
the two-forms 6" can be effectively assumed to be self-dual. As such they can be expanded in 
the background self-dual two- forms After such an ansatz is substituted into the action 

(I146p the result simplifies considerably. However, in order to exhibit the physical modes we 
need to introduce some convenient gauge-fixing. Inspecting (I16ip we see that it is possible to 
set to zero the following components of the b aa : 




(161) 



0. 



.6 



0. 



b 



0. 



(162) 
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This gauge turns out to be very convenient. We now write the gauge-fixed two- forms 6™ as 
follows: 



where the independent fields are now and 6° and the "strange" normalization 

coefficients are chosen in order for the Lagrangian to be obtained to have the canonical form. 

Substituting (11631) into (11461) and using the algebra of E-matrices we get the following simple 
effective low-energy action: 

gHigg. = _ J Q» b 5 + Q^_ + d^d^ + 0^0^ + 8^%^. (164) 

This form of the Lagrangian makes the reality conditions necessary to get a real theory 
obvious. Indeed, it is clear that the reality conditions are: 

(bir = bt (biy = bi (biy = bj, (&d* = & (m) 

These conditions can be compactly stated by introducing the following su(2) ®g valued object: 

b := (b 4 _ T + + b%) ®V++ (b%T_ + b 5 z T z ) ® (166) 
+ (b%T_ + b%) ®W + + (b 7 _T + + b%) ® W. 

and requiring it to be hermitian: 

b f = b. (167) 

The action can also be written quite compactly in terms of b. Indeed, using the pairing given 
by the (-, ■) metric in the Lie algebra we get: 

Cf f f s = -(d^,d,b) (168) 

for the low-energy E 2 k effective "Higgs" sector Lagrangian. It is thus clear that, at least 
in the low energy regime, the "Higgs" sector of our theory consists just of 4 complex massless 
scalar fields with the usual Lagrangian. It is not hard to show that in the finite k limit the 
content of this sector does not change, and is still given by massless fields. 
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7.7 Hamiltonian formulation for the Higgs sector 



In this subsection we obtain the Hamiltonian formulation of the sector spanned by V±,W±. 
After the analysis performed in the previous subsection such an analysis is not really necessary 
as we know what the propagating DOF described by this sector are like, and we even know the 
correct reality conditions. However, we decided to perform such an analysis for completeness, 
and also to confirm the reality conditions found from the Hamiltonian perspective. One finds 
the Hamiltonian analysis to be exactly parallel to that in the gravitational case, with even the 
final expression for the Hamiltonian being analogous. This subsection is quite technical and the 
reader is advised to skip it on the first reading. As in the case of gravity, we start by performing 
the space plus time split of the kinetic BF part. 

BF-Part. From (11441) our Lagrangian density is: 

c m, 9 s = 2V2 {&>*)" Yfc, _ (^y^r^-j 

-\ ((db 4 ndb 5 r(2v,u + ay + (^r^ 7 )^ - ) . (169) 

Now, denoting the indices 4, 5, 6, 7 collectively by a, we have: 

(db a f = -d b t ab , {db a ) a = d t aa - 2e abc d b b* c , (170) 
where we have introduced the configurational variables 

t aa := e abc b a hc . (171) 



We do not need an expression for the expanded Lagrangian (11691) because a more compact 
expression in terms of the conjugate momenta will be obtained below. For now let us compute 
the momenta conjugate to the configurational variables t aa . It is sufficient to compute just one 
of the momenta to see the pattern. We have: 

:= 9 -§j^ = ~ ^ (K b (dot 7b - 2e^d e bl f ) + Z 0a d b t 7b ) (172) 

-\{d,tl - 2e^d e bl f ) - \ {K b (dot 5b - 2e b ^d e bl f ) + Z z 0a d b t 5b ) . 

Comparing it to (I140p we see that 7r 4a = 2ia^. This is precisely analogous to the relation (11041) 
we had in the case of gravity. Indeed, the above relation can be rewritten as ir^ = — 2ig4 Q a", 
which generalizes (11041) . The other momenta are obtained as follows: 



-2ig aP al (173) 
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We now need to solve for the velocities in terms of the momenta and substitute the result 
into the Lagrangian. Similarly to the case of gravity the velocities can be obtained by taking 
the spatial component of the equations (11371) . We get: 

d tt - 2e^d b bl = iV2 E+ og + -L £+ % b + i £^ + \ E* b n 5b , 
d tl - 2e a bc d b b 5 0c = -iy/2 E - a 7 Q - -±= E~ \ %b - i E* a a* - ^ E* % , 

d tl - 2e a bc d b b 6 0c = iV2 E - a 4 + -±= E~ 6 7T 56 - i - ± E* 6 tt 7 6 , (174) 
dot 7 a - 2e a bc d b b 7 0c = -iV2E+ a a B - ±= E + \ Ab + i Sg a a^ + i E^ 6 vr 66 . 

The time projections of the equations (j!37j) are then the Gauss constraints. 

For the last step we start from a convenient expression for the Lagrangian. This is given by 
an analog of ( 11071) . which reads: 

Cff s = - 2g ab ^f b ay v (175) 



2V2E+ ^a^al + 2^2E~ ^ v a\a\ + 2E z ^aJa^ - 2E 2 " ,/ aJa£, 



where are the structure constants. Expanding it, and converting the spatial components of 
the connection into momenta we get: 

1 111 

-iV2Ej a (a^7r 7a - 7r 4a a^) + iv^Xo a (a^7r 6a - 7r 5a a£) 
+ i ^o a ( a o 7r 4a - vr 5a ao) - iSo a (ciovr 6a - 7r 7a ao) . 



We can now compute the Hamiltonian: 

WfiT = n aa d t aa - = 2n aa e abc d b bl + \g ab crg ps flpL hef 'K^'Ksf , (176) 

The obtained expression is not the full Hamiltonian. To obtain the later we need to add 4 
Gauss constraints that are obtained as the time components of the compatibility equations 
(I137p . We will not need an explicit form of the Gauss constraints since we already know from 
(I16ip what is generated by them. 

The BB-Part. Let us now consider the potential part (11591) . The corresponding Lagrangian 
density reads: 

L H B T = -KP-g^lfo*. (177) 
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Expanding the spacetime index we get: 

£ST = KQapWlL - \t a X) + iKgafiVSJ? 1 . (178) 

Total Hamiltonian. We now form the total Hamiltonian T-C Hl " s = H^p 9S — C^ 9S and 
integrate out the non-dynamical fields &Q a . We get the following expressions for these fields by 
solving their field equations: 

b a ,a = -9 a ^a bc d%- l -e a . (179) 

This should be compared with (11141) that we have in the gravitational sector. We now substitute 
this back to get the Hamiltonian with second-class constraints solved for: 

n m " s = lgab9 a V S fap^ bef V S f - i(e abc d b n ac )t a a + -g aP (e abc d b 7r ac )(e aef d e 7r Pf ) , (180) 

plus Gauss constraints with their corresponding Lagrange multipliers. Note also that the Hamil- 
tonian we have obtained is analogous to the one in the case of gravity (11161) . Indeed, there is 
similarly the ir 2 term and a (edii)t term with an imaginary unit in front. There is also a <9 2 7r 2 
term with a parameter of dimensions 1/M 2 as a coefficient. Note that for any value of the 
parameter k this Hamiltonian describes modes that are massless. To rewrite this Hamiltonian 
in terms of physical propagating modes we need to understand the gauge-fixing. 

Gauge-fixing. To choose a convenient gauge-fixing that eliminates the gauge transformation 
freedom, let us discuss what the two- forms £ , E 2 become after they get projected onto the 
spatial hypersurface. Thus, let us find analogs of relations ( 1971) . ( |96i) . Let us introduce the 
following 3 spatial vectors: 

£+ := m a , S 0a := m a , := n a . (181) 

Then, taking various projections of (11261) . it is easy to check that the following relations hold: 

m a m a = rh a rh a = 0, m a n a = rh a n a = 0, (182) 
m a fh a = 1, n a n a = 1. 

Taking different projections of (1 1 2 6 [) one finds the spatial pull-backs of the two-forms in terms 
of the vectors introduced: 

E+ b = n a m b - m a n b = -ie abc m c , £~ 6 = m a n b - n a rh b = -ie abc m c , (183) 

Y, z ab = m a m b - m a m b = -ie abc n c . 

We now use (I16ip to fix the gauge as in (I162p . In terms of the configurational variables t aa 
the gauge conditions read: 

t 4 - = o, t 5+ = 0, t 6+ = 0, t 7 ~ = 0, (184) 
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where our convention is that t a+ := m a t",t Q = m a t",t az = n a t"- 



Let us now find the consequences of the Gauss constraints. In terms of the introduced 
vectors m a ,rh a ,n a these read: 

d a t Aa - 4=m a 7r? - ~n a 7T 5 a = 0, d a t 5a + + = 0, (185) 

d a t 6a - ■j=fh a ir% + ^n a rf = 0, d a t 7a + ^=m a nl - ^n a < = 0, 

Introducing more compact notations 7r+ := m a TT^,TC~ = rh a 7i^, = n a ii^ and passing to the 
momentum space we have: 

i\k\& - -±=7T+ - = 0, i\k\t 5Z + i=7T 6 - + I** = 0, (186) 

i|A;|t fe - i=7r 5 - + = 0, i|fc|* 7 - + -^tt 4 + - 1*5 = 0. 

We now use these constraints to find the components of the momenta that are conjugate to the 
gauge-fixed variables (j!84p . We have: 

nt(k) = -iV2\k\t 7 * + ^=7T Z 6 , tt 5 -(A;) = iV2\k\t 6 * + -^tt*, (187) 

TTe (Ar) = -iv^t* - tt+(A;) = iV2>|t 42 - -^t£. 

Let us now substitute these expressions into the 7r 2 part of the Hamiltonian. Thus, we have for 
the first term in f 1 1 8 j) : 

- + 7rf(-*X(*)) - ^(-^(A;) + Trf (-*)<(*)) (188) 

+ ^«(-k)t 4Z {k) + 7l Z 7 (-k)t 7Z {k) - TC Z 5 {-k)t 5Z {k) - 7T Z 6 {-k)t 6Z {k)) 

-\k\ 2 {t 5z {-k)t 4z {k) + t 6z {-k)t 7z {k)) 

Let us now work out the second term in (11801) . We use: 

- i € abc = n a {m b m c - m b m c ) + m a {m b n c - n b m c ) + m a {n b m c - m b n c ) , (189) 
which can be easily derived from (11830 to write the second term in ( 1180[) as: 

i\k\(7r-(-k)r + (k) - n+(-k)t a -(k)). (190) 
Here we again passed to the momentum space and used 

d a (e ik X(k))=ik a e ikx t%(k), (191) 
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where k a = \k\n a is a vector in the direction of n a . This makes only two of the terms from 
(1 1 8 9 [) survive. Expanding and using the gauge-fixing conditions (j!84p we get for this term: 

i|fc|(7T4 {-k)t 4+ (k) + 7t?{-k)t 7+ {k) - 7T+(-k)t 5 -(k) - 7T+(-k)t 6 -(k)). (192) 

The total Hamiltonian in the E 2 -C K low energy limit is given by the sum of two terms, i.e., 
( DUD and ffT92]) . 

Reality conditions. Let us now discuss the reality conditions that are appropriate in the 
E 2 <C k low energy limit. It is clear that they can be determined by "completing the square", 
similar to what we have seen in the Hamiltonian formulation of the gravitational sector (in the 
low energy limit). Thus, let us write the total Hamiltonian as: 

n m ^ = -~ (vti-k) - ^^(-k)^ fa) - ^ 4z (ao) - ^\k\ 2 t^(-k)t^(k) (193) 

-~ (*?(-*) - ^ 6z (-^)) (*m - S^c*)) - \\k\ 2 n-k) t ^{k) 

_I (ttJ (-fc) -2i\k\t^(-k)) {n+(k)-2i\k\t i+ (k)) -2\k\ 2 t 5 -(-k)t 4+ (k) 
— (^7(- k ) ~ 2i\k\t 6 ~(-k)) (7r+(/fc) - 2i\k\t 7+ (k)) - 2\k\ 2 t & ~ (-k)t 1+ (k) . 

The form of the reality conditions is now obvious. Indeed, we introduce new momenta 
variables: 

2\\k\ 2\\k\ 
TfJ(fc) = K(k) + ^t 5z (k) } r 5 (k) = **(k) - -f±t**(k), (194) 

<{k) = - ^-t 7 '(k), Tff (fc) = ifik) + ^t^(k), 
Tf-(jfe) = ti-4 (*0 + 2i|A;|t 5 -(A;), n+{k) = n+(k) - 2\\k\t^{k), 
n+(k) = n+(k) - 2i\k\t 7+ (k), itr(k) = Try (k) + 2i\k\t 6 -(k), 

and then require the following reality conditions: 

TfJHO = -(Tfg(fc))*, *?(-*) = -(7fg(fc))*, (195) 

= -(*+(*))*, iff (-*) = -(*+(*))• 
t 52 (-A;) = -(t 42 (A;))*, t 62 (-A;) = -f/ 2 (A;))*, 
t 5 -(_A:) = _(t 4 +(A;))*, t e -(-k) = -(t 7+ (k))*. 

It is not hard to see that these conditions are the same as we have derived earlier in the 
Lagrangian framework, see (11651) . Indeed, the extra minus present in H195[) is due to the 
following transformation properties of the basic two-forms: 

(£+)* = -£- 6 , (£* 6 )* = -Kb (196) 
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that directly follow from (11831) . The obtained real positive definite Hamiltonian is that of 4 
complex massless scalar fields, so we have full agreement with our Lagrangian analysis above. 
Reality conditions and the Hamiltonian for the full finite k theory can be obtained via precisely 
the same method as in the gravitational sector case treated in the Appendix. We refrain from 
giving such an analysis in this work, as it becomes even more technical. 



7.8 YM sector 

In this subsection we work out the Lagrangian for the remaining part of the theory, which lives 
in the part of the gauge group that commutes with the background su(2). The total Lagrangian 
we start with is a sum of kinetic term ( 11471) and the potential term ( 11591) . with an extra sign 
in the potential term coming from the metric component gss = — 1- This gives: 

c ym = 2i e ^bl udp al + kP- ^b%b 8 pa . (197) 

The further analysis is greatly simplified by making use of the reality condition for the b 8 
two- form from the outset. Thus, as we will also confirm by the Hamiltonian analysis in the 
next subsection, the two-form b % needs to be purely imaginary: 

bl := -ib%, (1%)* = b%. (198) 

This immediately leads to simplifications as the real part of the Lagrangian (11971) is then given 
simply by: 

c ym = 2e^~bl v d p al - ^% v . (199) 
Taking a variation with respect to b we learn that: 

b% = -e^F", (200) 

where F pu = d p al — d v a^\ is the curvature of our U(l) gauge field, which is therefore, for real 
k, real. Substituting the result back into the Lagrangian we get: 

C YM = --{F pv f. (201) 

K 

This is the standard YM Lagrangian with the coupling constant: 

9ym = (202) 

To convert this into a physical coupling constant we recall that we need to multiply the La- 
grangian by 327rG, as this is exactly the prefactor that converts the canonically-normalized 
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graviton Lagrangian (19411 into the Einstein-Hilbert one. Thus, the physical coupling constant 
in our arising YM theory is given by: 

9ym = 4ttG«. (203) 

Realistic particle physics coupling constants are of the order of unity (and smaller), so we learn 
that the parameter n must be of the order M?, which is what we have been using in the previous 
subsections. 



7.9 Reality conditions for the YM sector 

In this subsection we perform the Hamiltoanian analysis of the YM sector with the main aim 
being to obtain the reality conditions used above. As in all other cases considered, the reality 
conditions become obvious once the Hamiltonian is written down. 

We start from the Lagrangian (11971) . Expanding 

e^blM = -2e ab % a d b al - t 8a (d a 8 a - d a a 8 ), (204) 
where t 8a := e abc bf c , we see that the momentum conjugate to the connection a 8 is 

n8a ■= i^nr = ~ 2it8a - ( 205 ) 

dd a s a 

The Hamiltonian is then: 

ri = 4ie b 0a d b a c - a d a vr + K(b b 0a + —n ttJ - -b 0a 7i . (206) 

We find the non-dynamical fields &Q a via their field equations and get: 

bl = --e abc d b a^+ 1 -rrl (207) 
K 4 

Substituting this back into (1206P we get the "physical" Hamiltonian 

n™ s = - k (e^d b a 8 c + |tt 8 °) ' + ^n 8 a . (208) 

It is now clear that the "correct" reality conditions that give rise to a real positive definite 
Hamiltonian is: 

lm(e abc d b a 8 c ) + ^Re(7r 8a ) = 0, Im(vr 8a ) = 0. (209) 



From (I207p and (12051) it is easy to see that these reality conditions are equivalent to the condition 

,8 



that b 8 two-form is purely imaginary, which is what we have used in the previous subsection. 
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Passing to the real phase space and imposing the Gauss constraint <9 a 7r 8a = as well as 
the transverse gauge condition <9 a a 8 = we get the following simple expression for the real 
Hamiltonian: 

n/YM ^ /o Hreal\2 , K /_ 8a\2 / 01 n \ 

™reaZ = -(^aOfe ) + — (tT ) , (210) 

which again confirms that the parameter k/8 plays the role of Qym- 



8 Interactions 

In this section we work out (some of the) cubic order interactions for our theory. Our main goal 
is to verify that the YM and Higgs sectors interacts with the gravitational field in the usual 
way, and that the YM-Higgs interaction is also standard. We start with general considerations 
on the cubic order expansion of our theory. 

8.1 General considerations 

The third variation of the BF-term is 

5 3 S BF = 4i J 3 5B 1 A [5 A, 5 A} 1 , (211) 
and the third variation of the BB-term is 

*S m = 4i f* x ( ^Jg^ (B SB)"(B SBfHB SBr» (212) 

+ 6 fl3?W B >"(« B,B '"l- (2i3 » 

As in the case of the quadratic order expansion, it is most laborious to compute the deriva- 
tives of the potential. We have already computed the second derivative above. The third 
derivative of V(h) is given by: 

d 3 V(h) _ 9IJ d 2 f | 9kl d 2 f | Qmn d 2 f 



Q h MN dh KL dh IJ n Q h MN dh KL n Qh MN dh IJ U dh KL 0h IJ 

Trh d 3 f 

dh KL dh KL dh IJ 
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where the third derivative of the function of the ratios is given by 



a 3 / 



dh MN dh KL dh IJ 
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(215) 



where stands for the derivative of f£ a with respect to its r argument and 



d 3 



TrhP 



P 



dh KL dh KL dh u \(Trh)P 



(Tr h)p dh MN dh KL 

,2 



V 



+ 



(Tr h)P +1 
p 2 (p+ 1) 



9u 



KL 
dh MN 



+ 9KL 



dh P u 



Q h MN 



+ 9mn 



dh P T 



1,1 



d~h KL 



(Tr h)P+ 2 
p(p + l)(p + 2) Trh p 



9u 9kl h p MN + gu 9mn h P RL + 9k l 9mn h v l3 x 



(Trh) 



p+3 



9u 9kl 9mn 



(216) 



The first derivative of a power of h IJ is given by (!6T|) . We have not found a sufficiently simple 
general expression for the second derivative of h^j 1 with respect to h MN h KL , but the expression 
( IBTT) can be easily differentiated for any given p. The above expressions can be used to obtain 
the third derivatives of the potential for our background. The results are given in the next 
subsection. 



8.2 Interactions with gravity 

In this paper we shall not consider gravitational sector self-interactions. They are easily com- 
putable, but since the main emphasis of this work is on unification, it is of much more interest 
to compute the interactions of other fields with gravity and their self-interactions. In this 
subsection we consider the coupling of non-gravitational fields to gravity. 
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Thus, at least one of the perturbation fields SB 1 is to be taken to lie in the gravitational 
sector. It is then easy to see that this is the only interaction in the cubic order. Indeed, 
where two of the three perturbation fields lie in the gravitational sector and there is only one 
non-gravitational perturbation, there is no interaction coming from the potential part since 



d 2 V(h) 



dh aa dh bc 



d 3 V(h) 



dh ea dh cd dh 



ab 



(217) 



where a stands for the non-gravitational part of the Lie algebra. There is also no interaction 
coming from the kinetic part of the action for the structure constant fj K is zero when two of 
the indices are in the su(2) part and only one index is in the non-gravitational part. Thus we 
only need to consider the interaction that is linear in the graviton perturbation. It is natural 
to expect that this coupling is that to the stress-energy tensor of our non-gravitational fields, 
and this will be confirmed below. 

The interaction coming from the kinetic term is only non-trivial for the Higgs sector fields 
(since the structure constant with two of its indices in the "YM" part of the Lie algebra and 
one in su(2) is zero since the YM and the gravitational parts commute). This interaction is of 
the schematic type h(db) 2 , which is as expected for scalar fields coupled to gravity. We are not 
going to work out this term, even though it is not hard to do it using the explicit formulas for 
the connections worked out above. 

Let us concentrate on the interactions coming from the potential part of the action, as being 
the most interesting one. The relevant derivatives of the potential are as follows: 



d 2 V{h) 



dh a Pdh ab 
d 2 V{h) 



(218) 



dh b ^dh aa 
d 3 V(h) 



z ~T-9ctP9ahi 



dh dp dh ca dh ab 



9a/3 



2(2i) 



Note that the fact that the first quantity is zero is not completely trivial, at it involves a precise 
cancellation of two otherwise non-zero terms. 



We can now compute the relevant interaction terms using (12 12ft . We need to divide this 
expression by 3! to remove the extra multiplicity introduced by taking the third variation of 
the action. An additional simplification comes from the fact that in the first term in the third 
derivative in (12 18[) we have a matrix projecting onto the tracefree part of the gravitational 
two-form perturbation matrix Sq^"^ . This part is zero when the parameter g — > oo, which is 
the limit of the usual GR that we are considering. Thus, this part drops out and we have for 
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the gravity-non-gravity interaction term coming from the potential: 




(219) 



Here the extra factors of 3 in the first term and 2 in the second come from expanding the 
general Lie algebra indices in (12121) . and the factors of i/2 come by using the self-duality of the 
background forms Xg^. To understand this expression it is useful to separate the coupling to 
the trace of the graviton perturbation, and to the tracefree part. Let us consider the trace first. 
Thus, we take: 

b% = ^oV> ( 220 ) 

with the field h being proportional to the trace of the metric perturbation h^. It is then 
easy to see that the expression ( 12191) vanishes on such gravitational perturbations. This is, 
of course, as expected, for both our YM and Higgs sectors are expected to be conformally- 
invariant (classically). Indeed, this is standard for the YM fields, and for the Higgs sector this 
expectation follows from the fact that the fields are (up to now) massless. Using (12 191) it is not 
hard to check that there is indeed no coupling to the trace part of the metric, which confirms 
our expectation. Note that this also provides quite a non-trivial check of our scheme, for the 
whole scheme would be invalidated if we had found that our YM fields couple to the trace of 
the metric. 

We now confirm that the coupling to the tracefree part of the metric perturbation is also as 
expected. We only need to consider the second term in (I219p . as the first term involves only the 
trace part of the metric perturbation. Let us consider the YM sector first. We now substitute: 

b% = Z a 0[ /h„ ]p , (221) 
and use the anti-self-duality of this two- form to get: 

dfL-YM = ~n»%vK pX Kbl a = -2^+^%^. (222) 

Here an extra minus is due to the metric on the Lie algebra. The physical Lagrangian is 
obtained from here by taking the real part. This makes only the term in the self-dual projector 
p+fiupx fasti contains the metric to survive. Substituting (12001) we get: 

41-™ = -e^F» v e^F af3 h Xa . (223) 
Expanding the product of two e's here we get: 

41-yM = -K p F ua h^, (224) 
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in which expression we recognize precisely the coupling to the stress-energy tensor that arises 
from the YM Lagrangian (12011) . The sign in front is different from that in f 1 2 1 j) because the 
variation of the metric with two upper indices is given by minus W . Thus, the arising coupling 
of our YM fields to the gravitational sector is correct. 

Let us now discuss the coupling of the Higgs sector to gravity. It is easy to see that in 
the low-energy approximation in which E 2 <C k and the two-forms b" are self-dual there is no 
coupling coming from the potential term. Indeed, we have already discussed that there is no 
coupling to the trace part of the metric perturbation. Thus, there is only the second term in 
(I219p that can contribute. However, it contains a factor of {^ pcrTX b b pr7 b^ x ) which is contraction of 
a self-dual Higgs two-form and an anti-self-dual gravitational one. So, it is zero and the only 
interaction term in the Higgs sector comes from the kinetic term of the action. As we have 
already discussed, it is of the h(db) 2 form, which is just the coupling of the metric perturbation 
to the stress-energy tensor of our set of massless fields. We are not going to work out the details 
as they are slightly messy, but we hope that the discussion given is sufficient to show that the 
interaction is as expected. 



8.3 Interactions in the non-gravitational sector 



Let us now concentrate on the interactions in the non-gravitational sector, most interestingly 
those between the YM and Higgs sectors. 

First, we note that there are no cubic interactions in the non-gravitational sector that come 
from the potential term. Indeed, such an interaction term involves three perturbation two-forms 
b" with the Lie algebra index outside of su(2). It is not hard to see that the corresponding 



derivatives of the potential vanish: 

d 2 V{h) 



dh^dh 



0. 



d 3 V(h) 



dh^dh^dh™ 



(225) 



Thus, at cubic order we only need to consider the interactions coming from the kinetic 
term. It is not hard to see that there are no self-interactions in the Higgs or YM sectors, but 
there are two possible types of interaction between these sectors. One of them comes from the 
term g a/3 b a f^a 1 a 8 , the other comes from b 8 f^aa a a^ , where a now stands for the Higgs sector 
index. The second of this is an interaction of the type (l/n)F(db) 2 , and is thus suppressed at 
low energies by E 2 /k. However, the first interaction is non-trivial and important even at low 
energies. In fact, it is not hard to show that this is the standard interaction of the gauge field 
a 8 with the conserved U(l) current of the Higgs sector that is charged under the YM subgroup. 
We are not going to spell out the details that are again slightly messy, but the important point 
is that the YM-Higgs sectors interaction is also as expected for a set of scalar fields charged 
under the YM gauge group (Higgs fields). 
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9 More general potentials: Mass generation 



Up to now we have for simplicity considered a very special class of potentials that depend only 
on the invariants constructed from the "internal" metric h IJ using the Killing-Cartan metric 
gij. It is not hard to show that due to the fact that the rank of h IJ is at most six, there is at 
most six such independent invariants, and thus only at most five ratios to be considered as the 
arguments of the function /(•) in (153|) . However, it is clear that these are not the only possible 
invariants. Indeed, the most general gauge-invariant function of h IJ can also involve invariants 
constructed using the structure constants fj K . For instance, let us consider 

ffhhh := f PQR f STU h PS h QT h RU , (226) 

where the indices on the structure constants are raised using the metric on the group. More 
generally, one can construct a matrix: 

(ffhh) IJ := f IQR f JTU h QT h RU (227) 

and build more complicated invariants from traces of powers of h IJ and (ffhh) IJ . This leads to 
a much more general set of gauge-invariant functions. In this section we shall study implications 
of such more general potentials. Our main point in this section is that these more general 
potential functions lead naturally to Higgs fields becoming massive. This is very important 
for phenomenology, for massless Higgs fields interacting with the "visible" YM sector in the 
standard way is obviously inconsistent with observations. 



9.1 Potential with an extra invariant 

For simplicity, in this paper we shall consider only one additional invariant given by ( 12261) . We 
shall see that such a potential is sufficient to generate masses for the Higgs sector particles. 
It is not hard to consider even more general potentials, but we refrain from doing it in this 
already lengthy paper. 

Thus, let us consider the potential depending on one more invariant 

V ih) = lLA F (lL>L,...,^im), (228) 
n \(Trh) 2 {Trh) n {Trhf) 

where we have divided (I226P by (Trh) 3 to make the potential homogeneous degree one. Then, 
the first derivative with respect to h is 

dm = 9u F + TrhdF_ 
dh IJ n n dh IJ ' 
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with (dF/dh IJ ) given by 

9F _y>, d ( TrhP\ F , d (ffhhh\ 
dh IJ p 8~h IJ \{Trh)p) n+1 dh IJ \{Trhf ) 

where F' is the derivative of F with respect to its argument (Tr h p / (Tr h) p ) , F' n+l is the 
derivative of F with respect to its last argument and 

d ( ffhhh \ 3 f PQ (I fj f s hp R h QS 3 ffhhh 



dh IJ \{Trhf) {Trhf (Tr h) 

Now, let us compute the second derivative of V with respect to h. We get 

d 2 V(h) g u dF g_KL_ dF_ Tr h d 2 F 
dh KL dh IJ ~ ~ dh KL + ~n~ ~&hF + ~n~ d~h KL d~h IJ ' 
with (d 2 F/dh KL dh IJ ) given by 



gu- (231) 



(232) 



d 2 F ^ , d 2 f Trh> \ d 2 (ffhhh) 

dh KL dh IJ p dh KL dh IJ \{Trh)P ) n+1 dh KL dh IJ \{Trhf ) 



TrhP 



dh IJ \ (Trh) 



f d ( Trh q \ d ( ffhhh\\ d 

+ V {WW) + Pin+1) dh^ \Jt7WJJ dh 1 

V ( F » 9 ( Tr ~ hP \ F» 9 ( ff ~ hhh \ \ JL_ I ff ~ hhh \ 

+ ^ \ {n+1)p dh« L {JtVW p ) + {n+1)in+1) dh^ \Jt7W J Jdh^ \(Trh)* J 

where F pq stands for the derivative of F' p with respect to its q argument and similar for F^ n+1 ^ 
and F'/ n+v ,f n+l y It is easy to show that 

Using the equations above, we obtain the following expressions 



dh ab 
dV 



d ~ h af3 

d 2 V 



=0, (235) 



K + Yj 9af3 , (236) 



dh^dh 

where we have set (F) = and defined 



K \ 

9ab 9a/3 + -^gcdfabfaP J ( 2 37) 



41 61 



(-^n.+lJO 



A = v n + iyu . (238) 
The parameter k is as before, see (I152p . with the function F(-) of one more invariant in place 



of/(0- 
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9.2 Higgs sector masses 



In this subsection we show that the new parameter A introduced above receives the interpreta- 
tion of mass squared of the Higgs sector scalar fields. To this end, let us work out the quadratic 
part of the action that comes from the potential, concentrating only on the A-dependent part. 
The ^-dependent part was already taken care of by setting the Higgs sector perturbation two- 
forms b" to be self-dual, and this is unchanged for our more general potential. Dividing fl5Tl) 
by two, using the self-duality of b°L in the second term and simplifying we get: 

Sf = ~ y J \gc*f^(K^)^tfJ - g«^% v . (239) 

We now substitute in this expression the expansions f)163p for our two-forms (in a specific 
gauge). It is not hard to see that only the term f^bfap contributes and we get: 

Sf = A J b 4 ~b 5+ + b 6+ b 7 ~ + b Az b 5z + b 6z b 7z = -m 2 Higgs {h\ b), (240) 

where 

m 2 Higgs = "A- (241) 

Thus, as all other physical parameters arising in our theory, the mass of the Higgs sector 
particles also comes from the defining potential. 



10 Discussion 



In view of the length of this paper it is probably appropriate to recap our logic and emphasize 
the main results that we have obtained. Thus, we have started with a generally- covariant 
gauge theory for a group G, with the action given by (pQ). At this stage all fields are complex 
and reality conditions are later imposed to select the physical, real sector of the theory. We 
then perform the Legendre transform and pass to the two-form field formulation ([6]). Our phase 
space analysis in section @] is only needed to get a better idea of what should be expected for the 
number of propagating DOF of the theory. It does not form an essential part of our argument. 
The main analysis starts in section E] where we analyze the simplest case G = SU(2) and show 
how it describes the usual gravity in the limit when a certain parameter of the potential is 
taken to be large, or, alternatively, for low energies. For a finite value of the parameter (or for 
Planckian energies) one gets a modified gravity theory with two propagating DOF. However, 
as the low-energy limit of our theory is still given by GR, we do not need to understand the 
nature of this modification for purposes of this paper. 

We start with the analysis of the pure gravity case by describing how the Minkowski space- 
time looks like in the language of two- forms, see ( l65l) . The action is then expanded to quadratic 
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order, and the field equations for the connection field are solved for, with the solution given by 
(1761) . After the solution is substituted into the action one gets the linearized kinetic term fl79l) 
as a functional of only the two-form perturbation. This is supplemented with the potential 
term part (IBBl . After the parameter g is taken to infinity one gets GR written in terms of 
two- forms, with a very compact linearized action (|79l) . This action is considerably simpler than 
the one in terms of the metric perturbation, and the relation between the two arises via (l90lh 
We also perform the Hamiltonian analysis of the linearized theory, to show how the usual two 
polarizations of the graviton arise in this language. In the g — > oo limit this analysis reproduces 
Ashtekar's Hamiltonian formulation of GR, in its linearized version. The main purpose of this 
analysis is to select the reality conditions for the gravitational sector. These are particularly 
clear in the Hamiltonian formulation, and later in the paper the same strategy of deducing the 
reality conditions from the form of the Hamiltonian is used for other fields. In this section we 
only discuss the rather simple reality conditions appropriate in the GR limit g — > oo. The finite 
g case reality conditions are deduced in the Appendix, for completeness. 

Once the SU(2) case is understood we enlarge the gauge group to G — SU(3). We take the 
same set of two-forms (1651) for the background, which thus selects in the sm(3) Lie algebra a 
preferred gravitational su(2) subalgebra. The analysis of the gravitational part is unchanged, 
but we have carried it out once more using a different basis in the Lie algebra (root basis), 
in preparation for the analysis of the non-gravitational sectors. These split into a part that 
commutes with su(2) and that will later be identified with the YM sector, and the part that 
does not commute with su(2) and becomes the Higgs sector. 

Let us start with the Higgs sector. As in the case of gravity, we first solve the equations for 
the connections a° in terms of the perturbation two-forms b° v and then substitute the result 
back into the action. The resulting kinetic part of the action as a functional of the two-forms 6" 
is given by (11461) . There is also the part (11591) coming from the potential. Similarly to the case 
of gravity, the role of the potential part, in the low-energy limit, is to set certain components 
of the two-form field 6" to zero. After this is done, the perturbation two-forms b® u becomes 
self-dual and can be expanded in the basis of self-dual two-forms ( l65i) . The coefficients in this 
expansion become our Higgs fields. They can be seen to be charged under the gravitational 
SU(2) subgroup, comprising two irreducible representations of spin 1/2 of SU(2). They also 
transform non-trivially under the part of the gauge group that does not commute with SU(2), 
and so they are not all physical. A convenient gauge is given by (11631) . Finally, our Higgs 
fields are charged under the part of the gauge group that commutes with SU(2), i.e. under 
the YM subgroup, which in the case of G = SU(3) is U(l). After a gauge is fixed, one 
obtains a Lagrangian for the physical fields, and this is found to be just the usual one for 
a set of 8 massless fields. We then determined the reality conditions needed to make it into 
a real Lagrangian with positive definite Hamiltonian. These can be read off from either the 
Lagrangian we have obtained, or from the Hamiltonian formulation that is also developed. The 
end result is a set of 4 complex (and at this stage massless) scalar fields with the usual real 
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Lagrangian (11681) . These fields are later made massive by considering a slightly more general 
set of defining potentials. 

We then analyze the YM sector, both in the Lagrangian and Hamiltonian frameworks. As 
usual in this paper, the Hamiltonian framework considerations are most useful for determining 
the reality conditions that need to be imposed. After these are deduced, the derivation of the 
Lagrangian becomes straightforward, with the result given by (I20ip . The YM coupling constant 
arises as (I203p . with the parameter k related to the first derivatives of the potential function 
via (11521) . 

We then discuss (cubic) interactions between the various sectors of our theory and confirm 
that they are as expected for such fields. Namely, the interactions of all fields with gravity are 
via their stress-energy tensor, and interactions of the Higgs sector with the YM fields is via the 
Higgs conserved current. 

Finally, we consider potentials more general than has been the case before, and show how 
the first derivative H238j) of the potential with respect to the new invariant becomes (minus) 
the mass squared (12411) of the Higgs sector fields. The parameter A = —m 2 Higgs can be both 
positive and negative so we have the possibility of the Higgs potential pointing both up and 
down, depending on the form of the defining potential. For negative m 2 Higgs and thus positive 
A the configuration b = is unstable and a new vacuum to expand about should be chosen, as 
in the standard Higgs mechanism. This finishes our demonstration of the fact that the content 
of the theory expanded around the Minkowski spacetime background is as desired. 

Let us now discuss whether the unification scheme described in this paper can be deemed 
"natural" in the sense that it naturally produces "realistic" values of the parameters such as 
masses and coupling constants. To this end let us see what dimensionful parameters are present 
in our theory. When the action is written in the form ([1]) the integrand has the mass dimension 
4 (assuming that the connection has the mass dimension one), and there are no dimensionful 
parameters in the theory at all. After the Legendre transform to ([6]) the two-form field has the 
mass dimension 2, and there are still no dimensionful parameters. However, since a part of this 
field is to be interpreted as the spacetime metric, it needs to be made dimensionless, and this 
is when a dimensionful parameter is introduced into the story. Rescaling the two-form field to 
give it the mass dimension zero introduces a parameter of the mass dimension 2 in front of the 
action (interpreted as 1/G, where G is the Newton's constant), as well as makes the potential 
function to have the mass dimension 2. This introduces a length (or mass) scale into the theory, 
and it is clear that there is only one natural mass scale given by M p . 

Various parameters of the theory are then obtained as derivatives of the potential function 
evaluated at the background, and these have mass dimension 2, or, in the case of the YM 
coupling (12031) as the product of the derivative of the potential times G. It is thus clear that 
the natural values for mass parameters arising in our theory are M p , and for the dimensionless 
parameters such as the coupling constant gyu ~ 1- However, this is precisely the values that 
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are "realistic". Indeed, as our Higgs fields interact with the "visible" YM sector, we need to 
explain why they are not observed. This is explained by their very high mass that makes 
them essentially irrelevant for the low energy physics. Second, the realistic values of the YM 
coupling constants of particle physics are order one, and precisely such values are natural in our 
unification scheme. Overall, our unification model is "realistic" in the sense that it reproduces 
everything that could be desired from such a simple setup. 

An important ingredient that is missing from our simple-minded model is that of the usual 
symmetry breaking mechanism of particle physics. Such a breaking, if present, would introduce 
additional mass scales into the theory and make it much richer. The model considered in this 
paper in which the bakground only broke the G symmetry down to the gravitational and YM 
ones did not break the YM gauge group. However, it is clear that our model naturally allows 
for such further breaking of symmetry. Indeed, we could take the background to be more non- 
trivial and give to some of our Higgs fields a non-trivial vacuum expectation value. Since our 
Higgs fields interact with the YM sector in the standard way, the effect of such a non-trivial 
VEV is also going to be standard - the YM symmetry is going to be broken, with some of the 
gauge fields becoming massive. It is then very interesting that in our scheme this standard 
particle physics symmetry breaking mechanism receives a new interpretation. Indeed, a non- 
trivial VEV for the Higgs is now on the same footing as a non-zero value for the metric. In 
other words, in our unification scenario the Higgs fields and the metric are just different parts 
of a single two- form field multiplet B 1 . Details of, for example, Hamiltonian analysis of the 
gravitational and Higgs sectors also confirm a very close analogy between the two. Thus, in 
a sense, it is the Higgs fields and the metric that become truly unified in our scenario. It is 
of considerable interest to study such more involved symmetry breaking scenarios. The goal 
would be to see if a truly realistic unification that puts together some GUT gauge group, a set 
of Higgs fields required to break it to the gauge group of the standard model as well as gravity 
is possible. This question is, however, beyond the scope of this paper. 

Yet another very important ingredient that is missing from our scenario is fermions. These 
are usually unproblematic for any scenario that operates in Minkowski spacetime. However, 
we start with a generally covariant theory with no metric in it, so it is not at all clear how and 
if fermions can be added. At the moment, this is probably the most serious objection against 
our scenario, but we remain hopeful that fermions can be described in our framework. The 
only possibility for this seems to be to further enlarge the connection field in ([1]) and make it 
" fermionic" . This might also require a " generalized" connection that is no longer a one-form, 
as fermions that we would like to obtain are not forms. We leave investigation of all these 
difficult but very interesting questions to further research. 

Finally, let us briefly touch on the question of quantization. The theory we have considered 
was classical, but, of course, it has to be quantized. It is then clear that our action §\§ is non- 
renormalizable in the usual sense of the word. Indeed, expanding the theory around Minkowski 
spacetime we have obtained a Lagrangian consisting of some renormalizable pieces - in the Higgs 
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and YM sectors - as well as gravity with its non-renormalizable interactions. However, there 
are also higher order interactions that are non-renormalizable, and the full action is given by 
an expansion containing an infinite number of non-renormalizable terms. Thus, the full theory 
is non-renormalizable. This is, of course, as expected, for we cannot hope to bring together a 
non-renormalizable theory (gravity) with renormalizable other interactions in a renormalizable 
unified theory. At best, we can hope for a non-renormalizable unified theory, and this is what 
is happening in our scenario. 

At the same time, what our starting action (CQ) describes is just the most general generally- 
covariant gauge theory. For this reason it can be expected that the class of theories obtained 
by considering all possible potentials /(■) is closed under renormalization. Indeed, all terms that 
could arise as counterterms are already included into ([I]) and so the only effect of renormalization 
should be in the renormalizing the defining function /(•). This expectation is consistent with 
the outcome of our analysis. Indeed, we have seen that, for instance, the YM theory gauge 
coupling is just a certain parameter of the potential defining the theory. This parameter is 
known to flow with energy, and from the perspective of our scheme this corresponds to a 
flow in the potential function. If the sole effect of renormalization is a flow in the space of 
potentials, the non-renormalizability of our theory ceases to be much of a problem. Indeed, 
it is then possible to hope for a non-problematic UV fixed point corresponding to some very 
special potential that would thus provide a UV completion of our theory. In this context it is 
interesting to remark that, since the gauge coupling is known to flow to zero value in the UV 
(asymptotic freedom), and such coupling in our scheme is on the same footing with e.g. the 
parameter g describing the strength of gravity modifications, it is natural to expect that g flows 
to zero in the UV as well. However, it is not hard to see that this corresponds to the defining 
potential V(-) flowing towards the one of the topological BF theory. Thus, at least prior to any 
concrete analysis, it seems that the sought UV completion may be given by the topological BF 
theory, something that in the past has been suggested in the literature in other contexts. All 
in all, the absence of the usual "finite number of counterterms" renormalizability of our theory 
may not be a problem as the theory may possibly be renormalizable in the sense of Weinberg 
[28J as containing all possible counterterms, see also [29j for a more modern exposition of the 
notion of "effective renormalizability". 

To summarize, there are many open problems of our scenario, notably questions of whether a 
realistic grand unification is possible, whether fermions can be described in the same framework, 
and whether the expectation of effective renormalizability is realized. However, it appears to 
us that in spite of all the open problems the scenario described already suggests some very 
interesting new interpretations and and is thus worthy of further exploration. 
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Appendix: Reality conditions for modified gravity 



The "correct" reality conditions for the full modified gravity theory can be worked out from 
the condition B l A (B^)* = 0. In linearized theory this becomes: 

S a A {b b Y = S b A b a , or S a ^(6^)* + S 6 ^fo^ = 0, (242) 

where (b a )* is the complex conjugate two-form perturbation and £ is given by fl67j) . We now 
rewrite this reality condition using the space plus time split. We get: 

i{t ab - {&)*) + 2{b ab + (6 fe a )*) = 0. (243) 

To get this condition we have used E£ c = e a bc , £q 6 = i5% and recalled the definition (|99|) of 
the configurational variable. We should now analyze this condition together with the already 
known solution (I115p . (11141) for the components b a Q. 

Let us first consider the trace and anti-symmetric parts of ( 12431) . Then in the tracefree 
symmetric gauge for t ab these conditions simply state that the lapse and shift functions N, N a 
are real. This explains why the factor of i was introduced in (11151) in front of the lapse. 

Consider now the symmetric tracefree part of (12431) . The corresponding components of 6 a g 
are known from (11141) and we arrive at the following condition on the phase space variables: 

^-Re (e^devr^ = lm(t a %. (244) 

In the case g — > oo that corresponds to GR this implies that {t ab )tf is real, but in the modified 
case the situation is more interesting. 

In addition to (12441) there is another condition that is obtained by requiring that (12441) is 
preserved under the evolution. Thus, we need to compute the Poisson bracket of ( 12441) with the 
Hamiltonian and impose the resulting condition as well. The computation is a bit technical, 
but at this phase space level there is no way to avoid it. Indeed, even in the case of GR it is 
clear from the form of the Hamiltonian (11161) that the relevant condition cannot be that the 
momentum is real, for the Hamiltonian would be complex due to the presence of the second 
term in the second line. The computation of the Poisson bracket can be done as follows. First, 
we introduce the real and imaginary parts of the phase space variables: 

j.ab j.ab i :j.ab „ab „ab < • ab fOAK\ 

t — t 1 + lt 2 , 7T = 7T 1 + 17T 2 . l^40J 

Second, we substitute this decomposition into the action written in the Hamiltonian form. The 
resulting action has real and imaginary parts. It is not hard to convince oneself that any one 
of these two parts can be used as an action for the system, the resulting equations are the 
same due to Riemann-Cauchy equations that follow from the fact that the original action was 



58 



holomorphic. We choose to work with the real part of the action. The relevant Poisson brackets 
are easily seen to be 

{nf(x),t lcd (y)} = 5i a 5 b j5 3 (x - y), (x)t 2cd (y)} = S^S^x - y), (246) 

with all the other ones being zero. The real part of the Hamiltonian (with the constraint part 
already imposed and dropped) reads: 

H real = h^ff _ hrf*> f _ e ^^\ f tf - e^TT^tf + — (Flftf - —{&"^f. 

We can now compute the Poisson bracket with the reality condition (12441) that becomes: 

-U^vrf = t f. (247) 
The Poisson bracket with the left-hand-side is: 

f H real efage bf} = __ A ab / 24g N 

The Poisson bracket with the right-hand-side is: 

{H real , tf} = ?rf + e efa d e t\ f - ~ g ^2 b - ( 24 9) 
Thus, the sought conditions that guarantees the consistency of (I247p is: 



+ e efa d e t\ } - ^Avrf = 0. (250) 



2g 

We now need to solve this for 7rf 6 , which gives: 

n f = _ e at i (251) 
2 l-A/2^' v ; 

where the denominator should be understood as a formal power series in powers of A/g. When 
j^oowe reproduce the GR result reviewed in the beginning of this subsection. 

We now have to substitute this, as well as the expression (I247P for into the action. This 
is a simple exercise with the result being: 

S real = J dtSx - \ ((tt&) 2 + (d a t b S R f)^J , (252) 



where we have defined: 



*GR = <, tGR = 1 _\ /2g - ( 253 ) 

These are the phase space variables in terms of which the Hamiltonian takes the standard GR 
form. This shows how an explicitly real formulation with a positive definite Hamiltonian can 
be obtained. We also see that for any finite value of g the graviton is unmodified. 
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